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ZUSAMMENFASSUNG. — Das im I. Teil dieser Arbeit zitierte WEINSTEINSChe Iterationsverfahren zur Strahldurchrechnung wird 


diskutiert. 


Summary. — The iteration method of ray-tracing described by Weinstein, cited in the first paper, is discussed. 


Sommasre. — L’auteur traite ici du procédé d’itération de WEINSTEIN, déja cité antérieurement, utilisé dans le calcul du trajet des 


rayons. 


4. Brechungs - und Ubergangsformeln. — Wir gehen 
aus von den im I. Teil dieser Arbeit abgedruckten, 


bekannten Brechungs-und Ubergangsformeln (5 a,b,c), 
und (6 a,b,c). Die Gleichungen (1) bis (7) des I. Teils 
stellen wir hier nochmal zusammen : 


(1) 2,=F,(q,), wobei g, = x + y% ist (Flachen- 


gleichung) ; 
; dew age 
( ) n= ix: dz ee 


(Scheitel-Kriimmungsra- 


dius) ; 


UST Aven tee 
Le le a Bey Gaz (s, ist die Sub- 
normale) 

eee = FT Ars Hy = Ht Ay; 
L ; 1 / dz(q) . 
feb) =C+As3;s = ~ 9/ dg’ | Brechungs- 

q=xwtty; formeln 
Ge): ¢' = yin’? Beye rae 2) 2 
ag, = 2, + 1,83 ys = Yr thm; 
(6b) 2, = 2, —@ + 4% 5 Ubergangs- 
(6c) 22 = F.(q2) 3 G2 = os ails Ys ) formeln 
(7) @ =A +l,—"N- 


* Berichtigung der Druckfehler im I. Teil dieser Arbeit : 


Im I. Teil dieser Arbeit (Opt. Acta, 1, 1954, Seite 127-140) 
sind folgende Druckfehler stehen geblieben : 


S. 129, GI. (3) : Statt « ¢ =... » muss es heissen : «G, =... » 
S. 130, zwischen Gl. (15) und Gl. (16) : Statt « ~7~ » muss 
@s heissen ; « —=~ ». 
z . 
S. 130, zweite Gl. (17) : Statt « (AzA — z,) » muss es heissen : 
# (Az — Az,) ». 
‘S. 130, Gl. (19 b): Statt «z= ... »m 
S. 133, rechts letzte Zeile : Statt «(z + Az)» muss es heissen : 
© + Az)» 
S. 139, unter dem Schema: Statt « fiir & z» muss es heissen: 
fir Tz». 


» muss es heissen: «Zz =... ». 
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2. Analytische Formulierung von Weinstein’s Ite- 
ration, soweit er sie geometrisch beschrieben hat. — 
Gemass der geometrischen Beschreibung, die WEIN- 
STEIN fur sein Iterationsverfahren gegeben hat, sucht 
man anstelle der direkten Auflésung der Ubergangs- 
formeln (6a, b, c) den Schnittpunkt der Geraden 
(6a, b) mit der Tangentialebene in dem Punkt D,_, 
der Flache (6c), der zuletzt als Naherungspunkt fir 
den Durchstosspunkt D ermittelt wurde, und man 
geht von diesem Schnittpunkt aus auf einer achsen- 
parallelen Geraden zu der Flache (6c) ber ; so kommt 
man zu einem neuen Punkt D, auf der Flache, der 
naher bei dem wahren Durchstosspunkt D hegt als 
der zuvor ermittelte Naherungspunkt D,—,. Mit D, 
als zuletzt ermitteltem Naherungspunkt kann man 
den im letzten Satz beschriebenen Prozess wieder- 
holen, um zu einem noch besseren Naherungspunt 
D+: zu kommen. 

Wir lassen in den Formeln (6a, b, c) iiberall den 
Index 2 weg, und wir bezeichnen mit (7; y; 2) die 
Koordinaten des zuletzt ermittelten Naherungspunk- 
tes fiir den Durchstosspunkt und mit (7; y; 2) die 
Koordinaten des neu zu ermittelnden Naherungs- 
punktes. 

Die Gleichung der Tangentialebene an der Flache 
(6c) im Punkt (a; y; 2) lautet : 


(8) (a — a) + yy —y) + s@—2) =0 ; 
dabei ist 

- 1 dz(q) 
a hs =|( dq eee pa 


Mit dieser Tangentialebene hat man die Strahl- 
Gerade (6a, b) : 


(10) pape; yHai tly 2=(4—el+le 


zum Schnitt zu bringen; die Schnittpunktskoordi- 
naten bezeichnen wir mit (x, ; ¥, ; 2,). Die beim voraus- 
gegangenen Iterationsschritt erhaltenen Schnittpunkts- 
koordinaten bezeichnen wir mit (2,3 ¥,3 Z.) : 
damit kann man die Strahl-Gerade auch so darstel- 
len : 


3-4 
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(11) caa,+ ke; y¥=y, thn; zz, + ke. 
Die vier Gleichungen (8) und (11) hat man nach den 
Koordinaten (2; ¥; 2) =(%,3 Yi 4.) des gesuchten 
Schnittpunktes und dem zugehérigen Parameter 
k =k, aufzulésen. Da man auf einer achsenparallelen 
Geraden von diesem Schnittpunkt (#, ; y, 5 4,) Zu dem 
Naherungspunkt (2; y; 2) fiir den Durchstosspunkt 
kommt, ist 

(12) « 


ebenso war &, aig Ve = y. Das beriicksichtigen wir 
bei der Auflésung der Gl. (8) und (14). Man erhialt : 


8(2 — 2,) 
(13a) a —; 
ze+ynt st 
(13b) e=a+tkhe; y=ythn; 
(13) Roney Pkyos 


Aus x und y erhilt man z und s durch die Flachenfunk- 
tion F : 


(13d) 2 = F@? + y’); 
(13 nth 
td =< 2 Fa? + 9) 

Die Gl. (13a, b, c, d, e) sind die Lterationsformeln fiir 
den Ubergang. 

Die Iteration beginnt damit, dass man die Strahl- 
Gerade (10) mit der Tangentialebene in irgendeinem 
zuerst angenommenen Naherungspunkt fiir den Durch- 
stosspunkt zum Schnitt bringt. Als ersten Naherungs- 
punkt kann man beispielsweise nehmen : 

A) den Punkt auf der Flache, der die gleichen late- 
ralen Koordinaten x und y hat wie der Durchstoss- 
punkt des Strahls an der vorhergehenden Flache 
(in den Ubergangsformeln (6 a, b) oder (40) mit dem 
Index 1 bezeichnet) — man legt dann dem ersten 
Iterationsschritt (43) die Werte. 


ty Ya Mae eae ee ee 


ae sel BS —4 
2=F@'+y*); s ~ 2a + 9) 
zugrunde ; 
oder 
B) den Flachenscheitel : z=0 -y =0;,=—r : 


man rechnet dann beim allerersten Iterationsschritt 
statt mit (13) mit den Formeln : 


(15a) L=— at 
(15b) 2=y+16s y=yntleq; 
(15) 2 =[4—e] + 46; 
(15d) & = F(a + y*); 
endl 
(15e) SO pi eianeeya Fi(a® + y?) 
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Nimmt man als ersten Niherungspunkt 
C) einen beliebigen Punkt der Fliche : 


ay 9; =F TY), 
—1 

2 FG + 9) 

man beim allerersten Iterationsschritt statt mit (13) 

mit den Formeln : 


mit der Subnormalen s = , dann rechnet 


a(%@—ay) + y(y—y1) + 5(2 —[2,—e]) 


(16a) 1, == = = 
ue + yn 4-55 
(16b) 2 =a +16; y=y tin; 
(16c) 2, =[%—e]+J,¢; 
(16d) 2 = F(x? + y%) ; 
= 
pi eerCn. 


Man konnte auch diese Formeln (16) bei allen Itera- 
tionsschritten anstelle von (13) verwenden; man 
braucht dann z, nicht zu berechnen und kann Gl. (16c) 
weglassen. Jedoch ist der Formelsatz (43) vorzu- 
ziehen, weil Gl. (13a) einfacher als Gl. (16a) ist, und 
weil man nach dem allerersten Iterationsschritt 
X13 Y¥13 [2,—e] bei den Formeln (43) nicht mehr 
bendtigt. — Die Formeln (16) gehen fiir az =10% 
y =0; 2 =—r in die Formeln (15) tiber. Ersetzt 
man in den Formeln (16), entsprechend dem Unter- 
schied zwischen den beiden Darstellungen (10) und (11) 
fiir die Strahl-Gerade, x, durch “3 y, durch y; 
[z, —e] durch z, und / durch &, dann erhalt man die 


Formeln (13). 4 
Ist die brechende Flache eine Kugelfldche, dann ist | 


(17') F(q) =— Vo 


folglich ist s = z, und aus den Iterationsformeln (13) 
wird : 


2(z —z,) 
(133%) k, = ——————_—¥ ;; 
rEt+yn+2e 
(13b’) e@=a+keé; y=ytho; 
(13¢’) 2, = 2, +k,6; 


(13d'‘) zg =— Vr? — (a + y®) | 


Im Fall A) legt man dem ersten Iterationsschritt 
(13’) die Werte zugrunde : 


(14°) 


eS US Ui Ba ee Oe 
é (2 + y? 


a r2 


— 
. 


(r) 


Im Fall B) rechnet man, ausgehend von z = 0; 


y =0; z=—r, beim allerersten Iterationsschritt 
mit 

F [2,—e]+r 
(15a’) bop : : 
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! Ae as 5 = . eI peal 

(15b") =a +les y=yit lo; a Man legt beliebige Werte &'; »' und dazu 

(15¢") Bee By G:| are aeas = + Vn’? — (B+ 7’) dem ersten Iterations- 


(15d’) oz = 


Vr 
(r) 


Im Fall C) rechnet man, ausgehend von beliebigen 


(eye). 


Werten 2; y und von z = aN pee (a? + 4?) , beim 
(r) 
allerersten Iterationsschritt mit : 
(16a’) l. a ae £4) + yly Y1) macs [ 2, é }) 
LE +yn ci 
(16b') ¢ =a,+16; y=ytlyn; 
M0e’) 2, =[4—e] +10 ; 
(16d’) z = —v/r* — (a? + y?) 
z () wos WEeos 


Wir haben noch die Iterationsformeln fiir die Bre- 
chung hinzuschreiben. Wir konnen sie ausden Itera- 
tionsformeln fiir den Ubergang zu einer Kugelfliche 
gewinnen, indem wir, entsprechend der Ahnlichkeit zwi- 
schen den Ubergangsformeln (6a), (6b), (6c) mit (47'), 
und den Brechungsformeln (5a), (5b), (5c), 
ersetzen : | 7 yy 21 ect al 
darchiis-2h7, 6 40 Bix! | +n! 


Wir erhalten folgende Jterationsformeln fiir die Bre- 
chung : 


wY 2 


(18) 


CYS | A 


c(t’ —%) 
(eG) ee 
B'a+ nyt t's 
Bop") 6 == Ene; 9) =! + xy; 
Mee’) 6, = %. + xs; 
(13d") a 35 a Mier (eS 7? n'?) ; 


Entsprechend dem Ubergang unterscheiden wir 
auch bei der Brechung drei Méglichkeiten, die Itera- 
tion zu beginnen : 

A) Man legt dem ersten Iterationsschritt (13”) die 
Werte zugrunde : 


(14") 


B) Man legt @’ =0; y' = 0; ¢’ =+ n’ dem ersten 
Iterationsschritt zugrunde, bei dem man dann statt 
mit (13”) mit den Formeln rechnet : 


15a") A= Pema 
$ 
‘15b”) 5 — é a I ae a! =» i mY 
(15c”) ee ae XS 
15d") v= —_ at Vn? vote (Gk ot. n'?) : 


(n‘) 


(n') 


schritt zugrunde, bei dem man dann statt mit (13”) 
mit den Formeln rechnet : 


(16a”) os ay: Be a 2) ae (x! 7) = AS = C) 
Ela aL n'y Be ele 

(16b") B= E+ ae; yi a=ant ray; 

(1607) C= C+ As ; 

(16d) [so ee. oe ep? (E'? + 9! 2) 


Nun diskutieren wir, wie schnell diese Iteration fiir 
den Ubergang oder fiir die Brechung konvergiert. Da 
fiir die Brechung gemass (18) das Entsprechende wie 
fir den Ubergang gilt, und da die Ubergangsfor- 
meln (13) und (15) nur Sonderfalle von (16) aa 
gentigt es, den Fall der Ubergangsformeln (16) z 
diskutionan’ Weichen x und y um kleine Gréssen i 
Ordnung » von den wahren Werten x und y ab, dann 
ist der Abstand des Punktes (x; y; z) von der Tan- 
gentialebene in (x; y ; 2) von (2)ter Ordnung klein. 
x und y weichen dann héchstens um kleine Gréssen 
der Ordnung (2 + 41) von den Koordinaten #, und Ys 
des Schnittpunkts des Strahls mit der Tampon. 
ebene ab. Weil z, = %; y, = y ist, haben also nach 
einem feceanonee tie die Naiherungswerte fiir die 
Lateralkoordinaten x, y des Durchstosspunktes nur 
noch einen Fehler der Ordnung. 


(19) = olds 


wenn sie vor dem Iterationsschritt einen Fehler der 
Ordnung » hatten. — Legt man dem ersten Iterations- 
schritt_entweder (A) die Werte «= 2,; y= He 
re = F(x? + y?) oder (B) die Werte x«=0;y¥=0; 

zZ = —r zugrunde, dann weichen x und y um Cresent 
der Ordnung # = 1 von den wahren Werten x und y 
ab, x und y also um Grossen der Ordung »’ = 3. Beim 
zweiten Iterationsschritt hat man o = 3; 0! =7. 
Beim dritten Iterationsschritt hat man o = 7; 
w' = 15. Beim N-ten Iterationsschritt hat man 


2% — 1) und. 
wo! = (Q8t?— 1). 


oa 
(20) 


Will man noch durch zusatzliche Naherungen die 
Wurzeln in der Iterationsformel (13d’) oder (16d’) 
fiir den Ubergang zu einer Kugelflache, und in der 
Iterationsformel (13d") oder (16d") fiir die Brechung 
wegbringen, dann hat man darauf zu achten, dass der 
Fehler, den man dadurch hereinbringt, die Konver- 
genz der Iteration nicht verschelchtert. Eine geeig- 
nete Niherung fiir die Quadratwurzel ist die NEw- 
Ton’sche Iterationsformel, die auch T. SmirH_und 
W. WEINSTEIN in ihren Formeln beniitzen : Ist w ein 


Niherungswert fiir Vo, dann ist 
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(21) 


1 (s v 
8a YN GS ene 
= EN W 


ein besserer Naherungswert fir Ve . Danach setzen 
wir fiir G]. (43d’) oder (16d’) : 

1/. r—(et+y') 

‘i ; 
Beim N-ten Iterationsschritt weichen x und y um 
Gréssen der Ordnung wo! = (28t!— 1) von x und y 
ab, 2 weicht daher nach Gl. (13d’) oder (16d’) von z ab 
um eine Grosse der Ordnung (#’ + 1) = 28+". Ent- 
sprechend weicht z von z ab um eine Grosse der Ord- 
nung (o + 1) = 2% = (m’ + 1)/2. Um den Fehler 
abzuschatzen, den die Niherung (22') gegentiber 
(13d’) oder (16d’) hereinbringt, setzen wir in die 
strenge Gleichung 


[c 
bo 
<7 
nN 
ES 
~ 
pales 


a 


rechts z = z 


2+ 
Ay Ye r2 —— (47 4 y?) 
Z= ta - -- = \= 
= 2 A 
aL 


Beim N-ten Iterationsschritt, wo z von z und damit 
auch von z um eine Grésse der Ordnung (@ + 1) 
abweicht, bringt also die Naherung (22') gegeniiber 
(13d’) oder (16d’) fiir z einen Fehler der Ordnung 
2-(@ + 1) =o' + 4 herein. So gross ist aber der Fehler 
von 2 gegentiber z schon nach Gl. (13d’) oder (16d’) 
wegen des Fehlers von x und y gegeniiber x und y, 
sodass wir, ohne die Konvergenz der Iteration zu 
verschlechtern, die Formel (13d‘) oder (16d’) mit der 
Wurzel, durch die Formel (22') ohne die Wurzel 
ersetzen koénnen. — Entsprechend kénnen wir auch 
die Forme! (13”) oder (16d”) ersetzen durch : 


lag ies at 
z | 


1s 
Dou i! = I 
aes oe [ 2 


3. Vergleich mit Weinstein’s Iterationsformeln. 
In unserer Schreibweise, ergiinzt durch die Bezeich- 
nungen 


(23) | eS Sa eM PNY Sie mae Teak, eM, 
lauten WEINSTEIN’s Iterationsformeln fiir den Uber- 
gang *) : 
Ce) OMIM Best irre wc: (Ps) 
plbatry)—35| + y°—2 p ds 
(W12) p, = a 


éx + ny + Cs(p,) 


1 Wir bezeichnen mit (W12) z. B. die bei WEINSTEIN, l. c. 
stehende Gleichung (12), umgeschrieben auf die Bezeichnungen 
der vorliegenden Arbeit. 
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(Ws) Los [ear Pilsner 
‘Sita ¢ | 
(oder 2, = [2, —e] + 1,8) ; ; 
| 


(W3) 2 =a tle; y=y tla; 


dabei ist die brechende Fliche dargestellt in der 
Form : 


(W1) s =s(p). | 


Wir vergleichen diese Formeln mit unseren Uber- 
gangsformeln (16). Die Gl. (W3) stimmen mit (16b) 
iiberein, ebenso Gl. (W4) mit (16c). Die Gleichnung 
(W12) stimmt aber mit (16a) nur naherungsweise | 
iiberein. Aus (16a) und (16c) erhaélt man durch Elimi- 
nation von /., und indem man beriicksichtigt, dass 


Ts y— z, [4 —@] H 
Rae i 1st : 
; z a . 
Ey 2) ts Z 
Z, = A acu i “ vis , und hieraus wegen (23) : | 
= ba + ny + Ss 
eo er) 
= Laas Nis gs 
(24) Heke ates vs = 4 
Ex + ny + bs 


Auf eine damit vergleichbare Form bringen wir jetzt 
WEINSTEIN’S Formel (W12). Nach (4) ist 


folglich 


4 P s(p) 
sr=—f sdp =—sp+ | pase 
2 0 s(0) 


Den zu p = p, gehorigen Wert von qg bezeichnen wir 


mit g*, er ist nicht gleich («? + y2). Damit wird 
aus (W412): 
eee (g; —49)+ 0/2 
ere — = ae = 5 
cS Sb + ny + ts(p,) 


Nun ist (q2 —yq) = 


=(48)--@.b+ 5 (54) - +00.) = 


(W12') 


=a ds = 
d S) a ee (+). s 0 s : 
und s(p,) = 8 pps ( Dg) OP Be) 
Damit wird aus (W12') : 
Fa (p —p,)-ts 
(W12") PSS SS = 
= bx + ny + bs 
C il ds 
4 + =—— =r SS = ss e 
bx + ny + ts 28 al 
-(p — p,) + O(p — p,)? 
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In der Abweichung der eckigen Klammer von 1 be- 
steht, wie man durch Vergleich mit Gl. (24) sieht, der 
Unterschied zwischen WEINSTEIN’s Iterationsformeln 
fiir den Ubergang und den unsrigen, die eine genaue 
analytische Beschreibung von WEINSTEIN’s geome- 
trisch beschriebenem Iterationsverfahren darstellen. 
_ Weicht nun bei einem N-ten  Iterationsschritt 
y von y (und & von x) um eine kleine Grosse a-ter 
Ordnung ab — wir schreiben dafiir x = x + ole) : 
y=y+ 0) —, dann gilt gemiss § 2 bei Anwen- 
dung unserer Iterationsformeln : 


B=at0 5 Fay p08); Zaz polo— 
(a5) 1 =c+0 5 y=y4+0) . 7 =z4 lot) . 

Cr LOCO, YY Oro ris Z,=2-+0%6) 

e =7+ 020+). y =y4Qot). ; = 74 e+e) , 


Wir bezeichnen den Wert von z, = p, —r, wie er 
nach Gl. (W12") herauskommt, mi z.(W), und die 
Werte von 2; y, wie sie nach Gl. (W12"), (W4), (W3) 
herauskommen, mit 2(W) ; y(W). Durch Vergleich 


der beiden Gl. (W12”) und (24) sehen wir : 
Ee tes 0W.6 (papa), 


Ss 


also wegen (25), wonach (p—p,) = (z—,) = glo) ist, 
(26a) ZW) 2 0h) 

Folglich ist nach. Gl. (W4), (W3) oder Gl. (16c), (16b) : 
(26b) 2(W) —2 = 0%); yw) —y = oh). 


Hieraus folgt wegen (25) : 


Bay a W) =a + 0h), y(W) = y + Ol, 


Nach WeErNsTEIN’s Iterationsformeln haben also 
nach einem Iterationsschritt die Naherungswerte fiir 
die Lateralkoordinaten 2,y des Durchstosspunktes 
einen Fehler der Ordnung 


(28) 


wenn sie vor dem Iterationsschritt einen Fehler der 
Ordnung » hatten. Betrachten wir den Fall, dass dem 
ersten Iterationsschritt solche Werte zugrundegelegt 
sind, dass und y von den wahren Werten « und y 
um Grossen der Ordnung o’ = 3 abweichen, dass man 
also eine paraxiale Naherung hat, dann hat man 
gemiss (28) nach dem zweiten Iterationsschritt 


wo’ =2e—1, 


’ =5, nach dem dritten Iterationsschritt o' = 9, . 


allgemein nach dem N-ten Iterationsschritt : 


(29) See ON Eye 


Mach WernsteINn’s geometrisch beschriebenem Itera- 
/onsverfahren und unseren daran anschliessenden 
}éerationsformeln ist die Konvergenz besser; es ist 
wach Gl. (19) : o =2+ 14, und nach Gl. (20) : 
ime = (271+ — 1). 

Zu dem gleichen Ergebnis kommt man auch bei 
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der Diskussion von WeinstTEIN’s Iterationsformeln 
fir die Brechung, die in unserer Schreibweise lauten : 


(W414) € = i. 
oy AG hate tee a) Ch Ere eae) 
(W5) a, = aaa 
Ss 


(W2) es = Ae n=ntay. 
Die Gl. (W5) und (W2) stimmen mit unseren Bre- 


chungsformeln (16c”) und (16b”) iiberein. Gl. (W14) 
weicht von (16a”) ab : Aus (16a”) un (16c”) folgt 


Stel ee 
Bat ny + Cs’ 
wahrend aus (W14) folgt : 


(30) 


(Wis) of —% = 


und hieraus : 

(W168) 
Coes Do one eae 

eo el (C.=5 +0 Ey] 
St+nyt t's 

Den nach Gl. (W14") herauskommenden Wert von 


bezeichnen wir mit ¢(W). Durch Vergleich der bei- 
den Gl. (W14") und (30) sehen wir : 


cw) — tf = ol). (@ — By, 


woraus entsprechend wie im Anschluss an Gl. (25) 
folgt, dass 


(31) e'(W) =e Q(20—1) . n/(W) ee Q2a—1) 


wird, wenn &! = #’ + 0!) 4% =! + 0) war. 


In der folgenden Tabelle ist nochmal die Ordnung 
der beim 1., 2., 3., 4., N-ten Iterationsschritt 
vernachlassigten Gréssen angegeben, vorausgesetzt 
dass der 1. Schritt eine paraxiale Naherung liefert ; 
und zwar bei den Iterationsformeln von T. Suirn, bei 
den Formeln von W. WEINSTEIN, und bei der geome- 
trisch beschriebenen Iteration von W. WEINSTEIN, 
fiir die wir in § 2 die Formeln angegeben haben. 


a! 2 3 4 N 
SMILE 2oyalaccnsisel cre 3 5 7, 9|j2N+1 
WEINSTEIN, Iterations- * 
Wormelin) scene eer 3 5 9117 |2>4+ 1 
WEINSTEIN, geometrisch reir 
beschriebene Iteration.|| 3 fi LS SSL 2: —1 
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Ebenso schnell konvergent wie WEINSTEIN’sS geo- 
metrisch beschriebene Iteration wire eine iterative 
Strahldurchrechnung mit den linearisierten Durch- 
rechnungsformeln des I. Teils dieser Arbeit ; diese sind 
aber als Iterationsformeln nicht so einfach wie die in 
diesem II. Teil angegebenen Formeln (13) und (13”), 
die WEINSTEIN’s geometrisch beschriebener Iteration 
entsprechen. Im J. Teil sind alle Brechungs-und Uber- 
gangsformeln linearisiert worden, nicht um zur itera- 
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tiven Strahldurchrechnung zu dienen, sondern zur 
analytischen Darstellung der Strahl-Aberration als 
Funktion der (LANGE’schen) Konstruktionsdaten 
und der Strahl-Anfangsdaten. Durch die Lineari- 
sierung ist eine solche Darstellung in Matrizenform 
moglich in einer Naherung, in der erst tertiére Bild- 
fehler (im monochromatischen Fall also erst Bildfeh- 
ler 7. Ordnung) vernachlassigt sind. 


Manuscrit regu le 15 janvier 1958. 


La diffusion résiduelle des surfaces polies et des réseaux 


A. MARECHAL 
Institut d’Optique, Paris 


SoMMAIRE. — On peut caractériser la qualité du poli d’une surface optique par Vimportance de ses déformations ainsi que par leur 
étendue. L’ étude de la lumitre diffusée peut fournir a ce sujet des renseignements utiles : on peut en effet exprimer 
— Vlénergie totale diffusée a UV aide du carré moyen des défauts, 
— la loi de répartition angulaire de U énergie 4 Vaide de la transformée de FourtER de la fonction d’autocorrélation des défauts. 
Le cas des défauts désordonnés de mise en place des trails des réseaux est également signaleé. 
Le présent article se limite aux développements théoriques, les vérifications expérimentales devant étre publiées aussit6t qu’elles 


seront completes. 


ZUSAMMENFASSUNG. 


- Man kann die Poliergiite einer optischen Fliche sowohl durch die Starke ihrer Abweichungen als auch durch 


deren Ausdehnung charakterisieren. Eine Untersuchung des Streulichtes kann in dieser Frage niitzliche Hinweise geben ; man 


kommt zu Ausdriicken fiir 


— die gesamte Streuenergie mit Hilfe des quadratischen Mittels der Abweichungen, 

— die Winkelverteilung der Energie mittels der FourtEeR transformation aus der Korrelalionsfunktion der Abweichungen. 
Auch der Fall der unregelmdssigen Fehler bei de Einstellung auf Gitterstriche wird in gleicher Weise behandelt. 

Der vorliegende Aufsatz beschrdnkt sich auf theoretischen Entwicklungen. Die experimentellen Nachpriifungen sollen mit- 


geteilt werden, sowie sie abgeschlossen sind. 


SumMARY. — The quality of polish of an optical surface can be characterised by the depth and extent of its deformations. A study of 


the diffusion of light can give useful information about this : 


one can express 


— the total diffused energy as a function of the mean square of the deformations, 
— the law of angular distribution of energy as a function of the Fourier transform of the auto correlation function of the defor- 


mations. 


The case of random-errors of position of grating lines is also treated. 
The present paper gives only a theoretical treatment ; experimental results in verification will be published when they are com- 


pleted. 


I. Généralités. — La diffusion intense des surfaces 
dépolies a déja été étudiée dans le cadre de loptique 
géométrique. Les irrégularités de ces surfaces grossieres 
sont en effet assez importantes pour que l’étude de la 
marche des rayons permette une interprétation accep- 
table des phénomeénes [1]. La discussion des expé- 
riences a, néanmoins, amené UrBANEK a la conclusion 
de la nécessité d’une théorie tenant compte de la 
diffraction, ce qui est indispensable lorsque les défauts 
de la surface deviennent petits. I] ne semble pas a 
notre connaissance qu’une telle étude ait été faite et 
il nous-a paru utile d’en présenter ici certains éléments. 

Cette étude concerne manifestement le cas des sur- 
faces optiques de bonne qualité, dont les irrégularités 
sont assez faibles pour que les déformations imprimées 
a la surface d’onde soient petites par rapport a la 
longueur d’onde. Nous montrerons d’ailleurs qu’un 
probléme analogue se pose dans le cas des réseaux, 
ou les déplacements accidentels des traits produisent 
une « diffusion » de l’énergie. 

Nous adopterons tout d’abord le cas schématique 
@une surface réfringente, présentant des défauts 


désordonnés mais de faible amplitude (fig. 1) et lon 
exprimera la déformation A de la surface d’onde &! 
en fonction de l’importance des défauts. On utilisera 


ae 


Fig. 1. 


a cet effet les techniques simples de l’optique géomé- 
trique. Si l’étendue des défauts n’est pas importante 
par rapport a la longueur d’onde un tel procédé peut 
paraitre criticable ; il serait facile en fait de montrer 
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qui est équivalent 4 l’application correcte du prin- 
cipe d’Huycens ov l’on prendrait la surface diffusante 
comme surface d’intégration, sans passer par linter- 
meédiaire de la surface d’onde. Remarquons d’ailleurs 
que si les irrégularités ont une étendue tres faible elles 
diffractent dans des directions trés éloignées de celle 
du rayon réguliérement réfracté (ou réfléchi) et il 
devient alors indispensable de tenir compte du carac- 
tere vectoriel des vibrations électromagnétiques. 
Notre étude se limitera en fait a la diffusion relative- 
ment proche. 


Il. Expression de A en fonction des défauts de la 
surface optique. Supposons pour commencer qu’il 
s’agisse d’une surface réfringente séparant deux milieux 
@indices n, n' (fig. 2) et soit I le point d’impact du 


Fig. 2. 


rayon lorsque la surface réfringente R a une forme 
correcte, I’ le point d’impact lorsque la surface, légére- 
ment déformée, est en R’. Le chemin optique L entre 
le point A de l’espace objet et le point A’ de l’espace 
image situé sur le rayon réfracté n’est pas le méme 
suivant que la surface est en R ou R’ et sa variation 
s’exprime par 
8h = s(n Al:-+ n' 1A’) . 
On peut exprimer commodément cette quantité a 
. . rat - . . 
Vaide des vecteurs unitaires u, w’ des rayons incident 
et réfracté 
>> > > > > > 
sh = nu 8] —n'u'.s1 = (nu—n'u') sl. 


D’aprés la loi de la réfraction, traduite sous forme 
vectorielle, on aura 
nu—n'u' = (ncosit—n'cosw')N. 
Si « est le déplacement du dioptre évalué sur la nor- 
a See : 
male N, on aura N.S! = « et il en résulte 


(1) sL = (ncosi—n’'cost’)e. 
La déformation A de la surface d’onde liée a L par 
la relation A = L, — L subira une variation 
8A = (n' cost’ —ncosl)e. 


Si en particulier le dernier dioptre d’un instrument 
est déplacé normalement a lui-méme d’une quantite ¢ 
-t dans le sens de propagation de la lumiere, on aura, 
puisque n’ = 1, 


8A = (cos i’ —ncosl)e. 
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Si maintenant l’on veut étudier le cas d’une surface 
réfléchissante immergée au besoin dans un milieu d’in- 
dice n, on aura, avec des notations évidentes sur la 
figure 3 (le déplacement de la surface étant cette fois 


Fie. 3. 


positif lorsqwil est effectué dans le sens de la lumiére 
réfléchie) 


> > > > > 
sL =n ds(Al + IA’) = n(u— uv’) 81 =22nN cos i SI 
ou encore 
(2) SA = 27 € COS 7, 


formule bien connue dans le domaine des interférences. 

Ainsi Vintroduction de légers défauts des surfaces 
réfringentes ou réfléchissantes produira des défor- 
mations irréguliéres de la surface d’onde : celle-ci ne 
sera plus sphérique mais oscillera de part et d’autre 
dune sphere moyenne, par rapport a laquelle elle 
présentera un écart fluctuant A. 


II]. L’énergie totale diffusée. — Nous imaginerons 
par exemple qu’une surface réfringente S approxima- 
tivement plane provoque dans la surface d’onde de 
petites déformations irréguliéres ; pour calculer tout 
d’abord Vénergie totale diffusée (en négligeant en fait 
celle qui serait diffusée sous de trés grands angles) 
nous imaginerons que !’on a placé une lentille L qui 
transforme la surface d’onde quasi plane en une 
surface d’onde quasi sphérique 2’ présentant les 
mémes irrégularités (fig. 4). Il est alors possible de 


iE 


SZ Lith 


Fic. 4. 


ealculer la répartition des amplitudes vibratoires en un 
point M(y’, 2’) du plan focal de la lentille L. On utilisera 
sur la surface d’unde 2’ les coordonnées angulaires 8’ 
y' qui sont les angles que fait la direction IF’ avec les 
plans w'F’z' et w'F’ y' (fig. 5). On sait que dans ces 
conditions on peut exprimer amplitude vibratoire 
(champ électrique) £(y'z’) au point M en utilisant une 
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Fig. 5. 


expression correcte du principe d’ HuyGens : si F(8’, y’) 
est la répartition des amplitudes sur une sphére centrée 


en F’, ona 
ai F(B!,y')e HB" 9! + 1'2" ag! dy! 


ou R = PF’ est le rayon de la sphére de référence, et k 
2m 


(3) Ely’,2') = 


est la constante k = 


Lorsque la surface d’onde est déformée d’une quan- 
tité A par rapport ala sphére de référence l’amplitude 
sur la sphére est de la forme Ey eWkA oy Ey est le 
champ électrique uniforme sur la surface d’onde. 
Puisque A est petit par rapport a A on peut écrire 
d’autre part 

(4) 

On peut encore choisir le rayon R de la sphére de 
référence de sorte que la valeur moyenne de A soit 
nulle : A = 0; dans ces conditions la répartitions des 
amplitudes sur la pupille est la somme de 

— une amplitude uniforme F, = Fy 


—une amplitude fluctuante et de valeur moyenne 
nulle 


(9) F,(8', x") = j k Ey A(B", y’) 

dans ces conditions le champ A (y', 2‘) est la somme de 
deux termes A = A, + A,; le champ A, est le champ 
diffracté par la pupille utilisée : il n’est appréciable 
que dans un domaine d’autant plus limité que la pupille 
est plus étendue; nous devons nous intéresser au 
contraire au champ A, dont la décroissance est beau- 
coup moins rapide : les petits défauts caractérisés par 
la fonction F, ont des étendues généralement petites 
par rapport a létendue de la pupille de sorte qu ils 
diffractent dans une région beaucoup plus étendue 
que celle de Ay. 

Pour calculer « l’énergie diffusée », i] faut tout 
d@abord en préciser la définition: nous adopterons a 
cet effet lintégrale du flux da a l’amplitude A, (bien 
que dans la petite région centrale ol A, n’est pas 
négligeable devant A, cette définition n’ait pas grand 
sens). 

En appliquant la relation qui exprime la conserva- 
tion de l’énergie on pourra écrire que cette énergie 
est proportionnelle a 


ff A,(y’, 2') |? dy’ da! =r ff F4(8", y') |? dB! dy’ 


et, en se reportant a la définition de /,, 


eA 1 4+ jka. 
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R? ke Bs {fs dp’ dy’ 


et Von fait ainsi apparaitre la valeur quadratique 
moyenne de A. Il est utile d’exprimer en fait l’énergie 
diffusée relative, en divisant par l’énergie totale égale 


ath Rs [foe dy’: 
ke I Ae dp! dy’ 
(6) — A 


jai? 


(fo 


Ainsi lénergie diffusée dans l’ensemble du plan 
image s’exprime aisément 4 l’aide du carré moyen 
de l’écart A (dans la mesure ot les défauts ont des 
dimensions suffisantes pour que la diffraction lointaine 
puisse étre négligée). 

ll est facile de relier ceci a la baisse d’éclairement au 
centre d’une tache de diffraction lorsqu’on introduit 
de faibles aberrations : on sait [2] que dans ce cas la 
perte relative d’éclairement est précisément égale a 
k2 A? et on comprend aisément que l’énergie ainsi 
soustraite Ala tache centrale se distribue dans le plan. 


IV. La distribution des éclairements. — Cette répar- 
tition est caractérisée par la fonction D,(y’, 2’) 
| A, (y’, 2’) |?etla fonction A,(y’, 2’) s’exprime (en vertu 
du principe d’Huycens) par la transformée de Fou- 
RIER de la fonction Fy. Or la fonction F, est aléatoire 
et il ne parait pas facile a priori de préciser autre chose 
a son sujet que la valeur moyenne de son carré. Si Von 
s’intéresse cependant a la répartition de l’énergie on 
peur remarquer que celle-ci est la transformée de 
Fourier de la fonction d’autocorrélation de F, [3]. 


Le principe d’Huycens se traduit en effet par la 


relation 


ik 1 a teal 
Olas Saul F,(8',y')e i889’ +2 )ap! dy! 


Larépartition de énergie s’exprime par | A,(y’,2')|? 
et on peut en exprimer la transformée de Fourier: 
si pw’ et v’ sont des fréquences spatiales dans le plan 
image, cette transformée de FourrEeR peut étre définie 
par la relation 


(8) d(u'.v') = I| A(y!,2') A¥(y’, 2! eRe 'y' +¥'2) day! dal, 
Cette fonction s’exprime aisément en utilisant un 
théoreme général sur l’analyse harmonique du pro- 


duit de deux fonctions et l’on obtient 


(9) d(u’,v’) = 
= re [fa y') F¥(g! J ' aot Ay de’ dy’. 


On arrive done a la conclusion suivante : une expeé- 
rience photométrique fournit directement la fonction 


| 


) 
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(y’, 2’) et si Yon effectue une transformation de Fou- 
RIER on peut connaitre d(u’, v’) qui est la fonction 
dautocorrélation de la distabuson aléatoire des F, 
(ou des A). 

Précisons maintenant comment on doit évaluer la 
répartition de lénergie diffusée lorsque la longueur 
d’onde varie: en écrivant F,(8', y’) = PEE A(p", ') 


ona 
GO) dia", v’) = 
eh? ki? Nie. B'—w'd, 7’ —v’a) dB’ dy’. 
Posons 
HA=B8B, wa=y et aes 


(qui est angle solide sous lequel on voit la pupille), 
on peut encore écrire la fonction d’autocorrélation de 
A sous la forme 


ay [fa (Sey) A(8 6, ¥ —y) dB’ dy’ =.SA7 Als.) 


ou A(8, y ) est une fonction d’autocorrélation norma- 
lisée : Sl ee = y = Oelle prend la valeur A(0,0) = 1; 
on a ainsi séparé les rdles de 

1°) amplitude des défauts, caractérisée par la 
valeur moyenne A? du carré de A, 

2°) ’étendue spatiale des défauts sur la surface 
diffusante : si les défauts ont une étendue notable la 
fonction A(®, y) est lentement décroissante : en effet 
Tes fonctions A(®6’, x’) et A(6’ — 8, y’—y) sont 
presque identiques et l’intégrale du produit est voisin 
de S A%. C’est le cas a) de la figure 6. Au contraire si 
A varie rapidement les deux fonctions A(p’ y’) et 
A(g’ — B, y’ — y) n’ont aucune relation et lintégrale 
de leur produit est nul : l’autocorrélation de A a une 
étendue trés limitée (cas 6 de la figure 6). La réparti- 
tion de l’énergie s’exprime finalement, en inversant la 
relation (8) qui définit la transformée de D? = | A,|? 


Paste = ff ae 


en tenant compte de (10) et (41), 
ee Rk? B28 M2 
[Adty's 21? = [faery 


Ba ee dy . 


On reconnait dans l’intégrale la transformée de Fou- 
RiER ®(Y.Z) de la fonction A, exprimée pour les 
valeurs Y = y'/Aet Z = 2'/a des variables ; on écrira 


Bale =) 


ee 
D,(y’,2') = 47? REGS tee, ieee 


e deme 'y’ +¥'2") du’ dy’ et, 


(12) 


(13) 


Ainsi la répartition de l’énergie diffusée ne dépend 


@ la longueur d’onde que par la combinaison 
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4 If 2! 

yee 

C NEN r 
car A n’en dépend que tres peu: la relation (1) permet 
d’en exprimer les variations en fonction de la disper- 
sion des divers milieux mais celle-ci reste généralement 
faible. 

L’étude expérimentale de la lumiére diffusée permet 
en principe de connaitre d’une part la valeur moyenne 
A? du carré de Aet d’autre part la fonction ®: on peut 
en effet mesurer pour plusieurs longueurs d’ondes le 


rapport entre l’éclairement D,(y’, 2‘) et le flux total 
F = R*S E?, qui traverse la pupille ; on peut ensuite 
4 


- ; x ; 
tracer la courbe qui représente = D,(y’, z') en fonction 


de y’/A et obtenir ainsi le produit A? ®; en fait la 
fonction ® est assujettie a la condition que sa trans- 
formée de FourteR A(8.y) soit égale a l’unité lorsque 

= y = 0 on devra done avoir, en posant y'/A = Y 
Clay Ave Z< 


I @(Y,Z)dYdZ =1, 


ce qui permet de connaitre séparément A? et ® a 
partir des données expérimentales. 

Pour fixer les idées a l’aide d’un exemple qui a des 
chances de représenter assez convenablement les 
phénomeénes, imaginons que la fonction d’autocorré- 
lation soit une fonction de Gauss (ce qui nécessite 
Vabandon de ’hypothtse A = 0, devenue inutile ici) 


Tae 
Ae 2 2 
) =e Bo 


A(B, 


ou By est un parametre qui fixe l’ échelle des corrélations : 
on peut dire que les irrégularités de la surface d’onde 
sont vues depuis le foyer F’ sous un angle de l’ordre de 
Bo. La fonction ®( Y, Z) s’écrit alors 


Lae Ze 
aces 


la décroissance de l’éclairement en fonction de y’ 
et 2’ s’exprime également par une fonction de Gauss. 
Il convient de remarquer que 

1°) siles défauts de poli de la surface diffusante sont 
de dimensions relativement importantes (fig. 6 @) G» est 
grand et la fonction décroit rapidement: l’énergie 
diffusée se localise pratiquement au voisinage immé- 
diat du faisceau régulier ; 

2°) si au contraire les défauts sont tres petits (leur 
corrélation est peu étendue) la lumiere diffractée 
s’étale largement dans le plan image (la fonction ® 
décroit lentement, ®) est petit), c'est le cas de la 
figure 6b. 

Cette remarque ne se limite d’ailleurs pas au cas 
particulier de la corrélation gaussienne mais a un 
caractere assez général : si la fonction d’autocorréla- 
tion A(®, y) s’étale dans un domaine important sa 
transformée de FouRIER sera au contraire rapidement 
décroissante. 
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Fia. 6a. 


Fia. 


V. Le cas du réseau. — Considérons un réseau plan 
comprenant N traits dont la mise en place peut étre 
affectée d’une erreur aléatoire, erreur commise sur la 
mise en place de chacun des traits étant indépendante 
de l’erreur commise pour les traits voisins. La somma- 
tion des amplitudes vibratoires pour une direction de 
diffraction quelconque s’exprime par 

N el P?. 

at 
Si l’on suppose que chaque trait, parfaitement bien 
mis en place diffracte une amplitude unitaire déphasée 
de l’angle » par rapport au précédent ; si maintenant 
la position du p™ trait est affectée d’une erreur e, ilen 
résulte une erreur de phase (pour un réseau par 
réflexion) 


19 c ‘ Fi 
$9, = Ep (sin Z =-'sin 2’) 


de sorte que amplitude sera, en fait, 


N 
», el Pei 5p (les déphasages 89, étant aléatoires). 
1 


On peut écrire amplitude, si les erreurs de mise en 
phase sont faibles, 


N 
Perv tite) =e re+ iY dress, 


1 


On sait que amplitude Ve! P ? est le plus souvent 
fond 


négligeable, sauf au voisinage des directions de 
diffraction réguliére. Par contre le 2 terme peut encore 


Avp,o) 


étre caractérisé par la moyenne de son carré (du fait 
que les phases p  prennent des valeurs trés variées) : 


J ape ee 
De Nee Ne 


et Véclairement obtenu est égal, en prenant pour 
unité léclairement central de la tache de diffraction, 


NES 
N2 


8 o., 


1 
AN 


expression qui peut conduire a des valeurs trés faibles 


a condition yue les erreurs sur les traits sotent effective-_ 


ment indépendantes (et ne présentent pas en particu- 
lier de périodicité, ce qui donne naissance a des ghosts). 

Signalons pour terminer que d’autres aspects de 
ces problémes ont été étudiés par A. LoHMANN et 
E. DrepERICHs avec qui nous avons eu le plaisir de 
discuter: 
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The interferometric measurement of metal oxide films 


J. > Kevey 
Physics Research Laboratory, University of Reading 


Summary. — An interferometric method is described for investigating oxide films in situ on the metal surface. The method yields the 
surface finish of the inside and outside surfaces of a metal oxide film, its volume ratio, its thickness and the location of the original 
surface with respect to the two oxide surfaces. All this information is contained in a single photograph provided the refractive 
index of the oxide is known. Some results for anodic oxide films on aluminitzn are given. 


SomMalRE. — Description d’une méthode interférométrique permettant d’étudier « in situ « des couches ad’ oxyde formées sur une surface 


métallique. 


La méthode donne V’état des surfaces interne et externe de V’ oxyde, son volume, son épaisseur et la position de la surface originale 
par rapport aux deux surfaces de l’oxyde. Ces données sont fournies par une simple photographie, si ’on connait déja Vindice de 


Voxyde. 


Quelques résultats sont donnés pour des couches d’oxyde anodique sur aluminium. 


ZUSAMMENFASSUNG. — Es wird ein interferometrisches Verfahren fiir die Untersuchung von Oxydschichten in ihrer Lage auf dem 
Metalluntergrund beschrieben. Das Verfahren liefert die Giite der Innen- und Aussenfldche der Oxydschicht, ihr Volumverhdltnis, 
ihre Dicke und die Lage der urspriinglichen Metalloberfldéche im Hinblick auf die beiden Oxydfldchen. Alle diese Angaben erhéilt 
man mit einer einzigen photographischen Aufnahme, wenn nur der Brechwert des Metalloxydes bekannt ist. Einige Ergebnnise 


fur die Anodenschicht aus Aluminiumoxyd werden mitgeteilt. 


I. Introduction. — When a metal is oxidised the 
oxide metal interface is below the original surface, 
step 5 in figure 1. If the volume ratio of the oxide, 
that is the ratio of the volume of oxide to the volume 
of metal from which it was formed, is greater than one, 
the top surface of the oxide will be above the original 
surface, step a in figure 1. 

Several workers have used inert markers to deter- 
mine the position of the original surface with respect 
to the two new oxide interfaces [ 1;2;3 ]. The technique 
is, however, usually applicable only to thick oxide 
layers or scales. It has the disadvantage of introdu- 
cing a third component into the metal, oxygen sys- 
tem. 

Other methods of film measurement which have 
been used measure only the total oxide thickness, 
a+b-+¢ (fig. 1). For thick films weighing may be 


REFERENCE FLAT 
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OXIDE 


used if the oxide density is known. Capacity measu- 
-vements have been extensively employed with the 

flm as dielectric ; here the dielectric constant and the 
-¢eue surface area must be known if the thickness is to 
/ be determined d. Surface factors of up to ten have 
| been reported, hence in any measurement ot thickness 
tne surface finish of the interfaces should be known. If 


the charge passed in forming an anodic film be used 
and a thickness calculated from Faraday’s law, we 
must know the efficiency of the anodic process and 
the oxide density. Other techniques for determining 
(a + b + t) include electron diffraction and anodic 
polarization [ 4]. 

The remaining methods are optical. Polarimetric 
methods in general involve laborious calculations 
although graphical methods may be used [5]. Fur- 
ther they require a knowledge of the optical constants 
of the uncoated metal. These are difficult to measure 
for metals which form, almost instantaneously, a thin 
protective thermal oxide film. The most simple and 
widely used method of estimating thickness is the 
comparison of the colour of the oxide with a standard 
colour gauge [6]. This standard gauge must first be 
calibrated by some other method such as weighing. 
This method is unsuitable for metals such as alumi- 
nium where the refractive index of the film is low 
and the reflectivity of the metal is high unless the 
oxide surface is coated to increase its reflectivity [ 7 ]. 
It requires a knowledge of the phase changes at the 
two oxide interfaces ; it is not accurate and like all 
the above methods does not give the location of the 
original surface with respect to the two new oxide sur- 
faces. 


II. Interferometric considerations. — Multiple beam 
interference has been used by Hass [8] to mea- 
sure the total thickness of a completely anodised 
evaporated aluminium wedge. It seemed to the author 
that much more information could be obtained by 
multiple beam methods. 

Since the layers to be studied are thin, FizEau 
fringes formed between the metal specimen and a 
reflector placed in contact with it were chosen. 
Although a greater accuracy may be attainable from 
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fringes of equal chromatic order [9], the dispersion 
of many of the oxides is not accurately known. 
Anodic oxidation produces the sharp steps which 
are most suitable for interferometric measurement, 
and has been used here. However, other methods of 
oxidation can be investigated by this technique. 
Figure 2 represents a specimen anodised over part of 


Fie. 2. — A: silver over thermal oxide ; B : silver 
over anodic oxide ; C : uncoated anodic oxide ; D: 
uncoated thermal oxide. 


its surface, areas B and C. Fringes running normal 
to the anodic boundary are obtained against a refe- 
rence flat. For a reasonably transparent oxide the 
fringes will be formed by reflections at the oxide- 
metal interface and the fringe step across DC will 
measure the oxide penetration, 5 in figure 1. If a strip 
across both the anodised and the unanodised surfaces 
is coated with a high reflectivity layer and fringes 
similarly obtained, areas A and B in figure 2, then the 
fringe step will measure the height @ in figure 1. 

Step a is measured in air and is obtained directly 
from the fringe displacement. Step 6 is measured in 
the oxide and hence its refractive index is required. 
The path difference for the step 0 is in figure 1. 


Aa[2={d+ (a+ b t t)n}—lat+d+nt} 


(A 2/2) —a(n—1) 
n 


or a 


where x is the fringe fraction measured and 4 the 
monochromatic wavelength used. It is assumed here 
that the index of the thermal oxide is the same as that 
of the anodic oxide. If this should prove not to be so 
allowance can easily be made for it. 

In the measurement of step a the reflecting surface 
is the same on both sides of the step and is a silver-air 
interface. Similarly, for b where the interface is oxide- 
metal. There is therefore no differential phase change 
provided the silver layer is thicker than 500 A [10] 
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and the structure of the surface is the same. By exa- 
mining the fringes formed on the four areas used for 
measurement in figure 2 we have a measure of the 
surface roughness of the top (A) and bottom (D) sur- 
faces of the thermal oxide, and the top (B) and bot- 
tom (C) surfaces of the anodic oxide. 

When comparing the surface roughness of the four 
areas it must be remembered that any roughness at 
the metal-metal oxide interfaces, C and D, is measured 
in the oxide and hence allowance for its refractive 
index must be made. 

The volume ratio of the oxide is clearly from 
figure 4 
a+b 


Q= 
b 


All this information can conveniently be recorded 
on a single photograph. To determine the thickness of 
the film (a + 6 + t), we need to know t which can be 
obtained from a subsidiary interferometric measure- 
ments, or in the case of some oxides, from the published 
measurement by other methods. It is normally 
small, t < 50 A. 


III. Experimental method. — If the spacing bet- 
ween the two interferometric surfaces is kept small 
large tolerances in collimation and line width of 
source are permissible. The reference flat was hence 
placed in contact with the silvered anodic specimen. 

The sharpness of a Fizeavu fringe system in reflec- 
tion is particularly dependent on low absorption [ 11 ]. 
Therefore a dielectric multilayer reference surface 
was made. This has the additional advantage that it 
is much less subject to attack by immersion liquids 
than a silver surface would be, in cases where it is 
desirable to fill the gap between the plates with a 
liquid of the same index as the oxide. Further, since 


the high reflectivity of the multilayer occurs only at a - 


specific wavelength, fringes may be obtained at this 
wavelength, and then the surface examined normally 
at a wavelength for which the multilayer reflectivity 
is low, [ 12]. 

The conventional type of multilayer made of suc- 
cessive quarter wave layers of magnesium fluoride 
and zine sulphide was not durable enough to stand 
the repeated application and adjustment of metal 
specimens. Ceric oxide proved to be a successful high 
index material when a suitable evaporation technique 
was developed. A seven layer multilayer of alternate 
quarter wave layers of ceric oxide and magnesium 
fluoride has been used for months and cleaned repea- 
tedly without any signs of damage. A Vickers Projec- 
tion Microscope was modified to enable multiple beam 
fringes to be obtained. The mercury green line, 4 = 
5461 A, from an under-run high pressure mercury 
lamp was used as source. 

The specimens were partially coated with a high 
reflecting layer. Silver evaporated through a series of 
bevelled slots was usually used. This gives a number 
of regions such as figure 2 on the one specimen. 
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The anodised specimen was placed on the reference 
flat and the fringe spacing and orientation adjusted 
by applying slight pressure to the specimen, or by 
placing a small weight on it. 


IV. Results. — Figure 3 shows the fringe pattern 
obtained with an anodic film grown on evaporated 
aluminium at 4.5 mA/cem? to 90 V in an ammonium 
borate electrolyte. The four regions are lettered as in 
figure 2. The reduced sharpness of the CD fringes is 
due to the lower reflectivity of aluminium compared 
with the silvered surface AB. The fine structure of 
the four surfaces is clearly visible in the fringe pattern. 
For this specimen a = 280 A, b = 884A andQ = 1.32. 


Fic. 3. — Anodic film grown to 90 V at 4.5 mA/cm? 
on aluminium. a = 280 A, b = 884 A, Q = 1.32. 


A film grown to 90 V at 14.3 mA/cm? is shown in 
figure 4. The high current leads to a more pronounced 
edge effect. For this specimen a = 294 A, b = 957 A 
ona Oo== 1.31. 

If care is not taken when anodising, the electro- 
lyte may be splashed or creep up on to the specimen 

leading to edge effects such as in figure 5. This can 
1 reduce the accuracy of measurement particularly if 
| tne specimen is slightly curved as this one is. This 
‘specimen was anodised to 60 V at 13.4 mA/cm? and 
[tae results were a = 172A, 6 = 488 A and Q = 1.35. 
- Typical results of measurements made of anodic 
exide on polycrystallic evapored aluminium films are 
@ ven in Table 1. The anodic layers were grown at the 
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Fic. 4. — Anodic film grown to 90 V at 14.3 mA/cem? 
on aluminium. a = 294 A, b = 957 A, Q = 1.31. 


constant current density shown in column 2 to a final 
voltage, the anodic voltage, shown in column 1. 

In the calculation of 6 and Q, Hennia’s[ 13 | value 
for the refractive index of anodic aluminium oxide 
has been used. For 5461 A this is 1.563 + 0.005. 

Each measurement in Table I is the mean of at least 
ten observations. The standard deviation for the fringe 
fractions a and 2 2/2 is normally less than 5 A. 

The volume ratio for anodic oxide grown in an 
ammonium borate electrolyte on evaporated poly- 
crystalline aluminium is seen from Table I to be inde- 


Fic. 5. — Anodic film grown to 60 V at 13.4 mA/cni2 
on aluminium. a = 172 A, b = 488 A, Q = 1.35, 
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TABLE 1 oxide interfaces with respect to the original surface. 
— = Panera _. The only additional information required is a know- 
Ie | Volume ledge of the refractive index of the oxide and, if total 
Nemes | Pats al? ue oe | pony : thickness is required, the thickness of thermal oxide 
Toltace eam) ~~ | 
; Voltage | mA/em? Q=—— | on the untreated metal. 
| [pei amet ae. 73K es A limitation is the need for a surface of good opti- 
39 18 |. 5371, “90 307 1.32 cal quality. This is not difficult to achieve on bulk 
| i specimens as with microinterferometric technique 
of see e029 ee si soe the flat area required is small. 
60 3.8 987°") 0-179 566! 13289 | | 5 
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Etude des angles de raccordement des surfaces liquides avec les solides 
par les interférences en lumiére polarisée 


M. Frangon et Mile Y, Gannon 


SommalRE. — Lorsquw un liquide est mis en contact avec un solide, 


par suite des forces d’adhérence et des forces de cohésion,la surface 


liquide se reléve ou s’abaisse suivant que le liquide mouille ou ne mouille pas le solide. L’angle formé au contact du solide par la 
surface liquide avec la surface solide est angle de raccordement. 


Nous nous sommes proposés au cours de ce travail d’appliquer une méthode d’interférence en lumiére polarisée & la mesure 


des angles de raccordement. Cette méthode, basée sur le dédoublement de V’onde incidente 


déterminer en chaque point, la pente de la surface d’onde. 


par un systéme biréfringent, permet de 


Nous nous sommes essentiellement attachés a définir les possibilités d’utilisation de cette méthode dans le cas de la mesure des 


angles de contact et a rechercher Vallure des phénoménes dans quelques cas particuliers : 


driques immergés. 


goutles de liquides, lames ou fils cylin- 


ZUSAMMENFASSUNG. — Bringt man eine Flitssigkeit mit einem festen Képer in Berithrung, so hebt oder senkt sich der Fliissigkeits- 
spiegel infolge der Adhdrenz- und Kohdsionskrdjte, je nachdem, ob die Flissigkeit den jesten Kérper benetzt oder nicht. Der Winkel, 
den der Flitssigkeitsspiegel gegenitber der Oberfldche des festen K6érpers bei der Berithrung bildet, ist der Berithrungswinkel. 


Wir schlagen nun in dieser Arbeit ein Interferenzverfahren im polarisierten Licht vor, 


um den Berithrungswinkel zu messen. 


Diese Methode beruht auf der Verdoppelung der einfallenden Welle durch ein doppelbrechendes System und erlaubt, die Neigung 


der Wellenfldche in jedem Punkt zu bestimmen. 


Wir haben uns in erster Linie bemitht, die Anwendungsmé 
winkel festzustellenund den Verlauf der Erscheinungen in einig 


Platten und zylindrischen Fédden. 


glichkeiten dieses Verfahrens fiir die Messung der Beriihrungs- 
en besonderen Fdallenzu untersuchen, z. B. an Flissigkeitstropfen, 
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SuMMARY. — When a liquid is placed in contact with a solid, the liquid surface rises or sinks according as to whether the liquid wets 
or does not wet the surface. The angle formed by the liquid surface at the contact with the solid surface is the angle of contact. 

We have applied an interference method using polarised light to the measurement of angles of contact. This method, based on 

the doubling of the incident wave by a birefringent system, makes possible a measurement of the slope of the wave surface at each 


point. 


The main purpose of our work was to determine the possibility of using the method for measuring angles of contact, and to 
explain the phenomena in particular cases : drops of liquid, cylindrical films or threads. 


Dans l’étude du mouillage d’un solide par un liquide, 
les chercheurs se proposaient autrefois de mesurer 
directement l’angle de raccordement. Etant données 
les difficultés rencontrées, les mesures de ce type, qui 
sont des mesures statiques, ont été remplacées par 
des mesures dynamiques parmi lesquelles nous cite- 
rons la méthode de GuAsTALLa. 

Dans la méthode nouvelle que nous proposons, la 
mesure de l’angle de raccordement peut se faire sim- 
plement méme dans des régions d’étendues extréme- 
ment faibles que lon peut observer au microscope. 
Les variations de l’angle de raccordement au contact 
d’un solide, variations dues a des causes diverses, 
peuvent étre étudiées et mesurées point par point 
avec précision. C’est en cela que la méthode pourra 
étre utile et donner des informations plus précises sur 
le mécanisme du mouillage. Elle permet de détermi- 
ner non seulement l’angle de raccordement mais 
également la tension superficielle. Quoique la mesure 
des valeurs absolues de ces grandeurs caractéristiques 
du mouillage soit facile, la méthode nous semble inté- 
ressante surtout dans |’étude et la mesure des varia- 
tions, locales ou non, de ces grandeurs. 


I. Principe de la méthode employée. — On se pro- 
pose d’étudier la déformation d’une surface liquide au 
contact d’un solide. Le liquide est en O dans la cuve G 
(fig. 1) et on étudie sa déformation au voisinage de la 
lame solide P immergée. Le systeme étant éclairé en 
lumiére paralléle, onde incidente est l’onde plane 2}. 
Aprés traversée du liquide la surface d’onde est défor- 
mée et devient l’onde Z,. La lumiére traverse ensuite 
un polariseur & et un systéme biréfringent I. qui 
dédouble l’onde X, en deux ondes ©, et 4. Ces deux 
ondes sont décalées une par rapport a l’autre dans 

deux directions : 

a) dans la direction de propagation de la lumieére par 
suite de la différence de marche introduite a la traver- 
-sée de la lame biréfringente. 

b) dans une direction perpendiculaire par suite du 


(1] (2) 


| 
| 
| 
Lo | 
| 
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dédcublement des ondes produit par la biréfringence 
du systéme. 

Les vibrations correspondant aux deux ondes X, 
et X, sont cohérentes par suite de la présence du pola- 
riseur &. Elles deviennent paralléles grace a lanaly- 
seur ,\, placé apres le systeme biréfringent L. Les deux 
ondes , et Z, sont reproduites a plus grande échelle 
sur la figure 2. 
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Supposons que dans les régions (1) et (4) ot ces ondes 
sont paralléles, leur difference de marche 8 soit égale 
a une longueur d’onde (jaune du sodium). Si on opere 
entre polariseurs croisés, on observe dans les régions (1) 
et (4) le pourpre du premier ordre. Par contre dans les 
régions (2) et (3) par suite de la déformation de la sur- 
face d’onde et du dédoublement latéral produit par 
le systeéme biréfringent L, la différence de marche 
nest plus égale a 8. Dans les régions (2) et (3) les 
teintes ne sont pas les mémes que dans les régions (1) 
et (4) ot les ondes sont planes. Si le décalage latéral 
n’est pas trop grand, l’ensemble des deux ondes &, et 
x, moule lobjet avec une précision suffisante pour 
que Vimage donne une représentation conforme de 
Vobjet. On choisira un décalage latéral de ordre de 
grandeur du pouvoir séparateur du systeme d’obser- 
vation. Dans ces conditions, il n’y aura pratiquemeut 
pas de différences entre ensemble des deux ondes ©, 
et , et la surface d’onde 2, de la figure 1. La méthode 
permet donc de déceler par des variations de teintes 
les déformations de la surface d’onde %, c’est-a-dire 
les déformations de la surface liquide. 


II. Détermination de la pente de la surface d’onde 
en chaque point. — Considérons deux portions assez 
petites des surfaces d’onde %, et Z, pour pouvoir 
étre confondues avec leur tangente (fig. 3). d représente 
le décalage latéral de ces deux ondes. La différence de 
marche PO = A entre les deux ondes peut s’exprimer 
en fonction de d et de la différence de marche 8 entre 


eae 
PLL 


les deux ondes X, et &, dans les régions ow elles sont 
planes. Ona: 


(1) - A= PO-= PN NO) =oeoedee 


L’examen de la teinte détermine par conséquent la 
pente « de la surface d’onde. En chaque point la 
teinte permet d’évaluer A. D’autre part, on peut 
observer la teinte dans les régions ot: onde est plane 
(pas d’objet). On a alors de la méme fagon 8. Connais- 
sant A et 8 on a la pente « correspondant a A. 


x 


III. Application 4 la mesure des faibles angles de 
raccordement. — La figure 4a montre une goutte 
de liquide sur une surface plane horizontale et la 
figure 4b représente le ménisque au contact d’une 
paroi solide verticale. Considérons par exemple la 
figure 4a. L’angle 8 est angle de la tangente a la sur- 
face du liquide avec Vhorizontale A’x. Soit le rayon (4) 
qui traverse le ménisque en un point quelconque et le 
rayon (2) qui le traverse dans la région ot |’épaisseur 


em 
dice du liquide, le chemin optique entre A et A‘ est 


Ay = Ne 4 6 CN re eee 
On a d’autre part, puisque 8 est supposé petit, 
de 


Or 


Fig. 4, 


est maximum. Si A’P est égale a e, et sin est Vin- 
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D’ot 
dA 
(2) a = (Ne 1G, 
dx 
dA, 


a represente la pente de la surface d’onde au point 
P étudié. Si « est cette pente, on aura 
(3) (t=), B: =e 


On a vu comment on pouvait déterminer « d’apreés (1). 
On aura donc 


4 ee ee 
(4) B aya 


En mesurant la différence A — 8 tout A fait au voi- 
sinage de M, on aura Pangle au contact cherché. 

Dans le cas d’une lame solide L immergée et incli- 
née (fig. 5), @ n’est pas égal a langle de contact et si J 


représente l’angle d’inclinaison par rapport a la verti- 
eale, on aura 


(5) 6 i eegi er 16 


La formule (4) permet donc de déterminer langle 
de raccordement 9. I] faut pour cela mesurer la diffé- 
rence de marche A — 8 correspondant a la région M 
des figures 4 et 5. Cette mesure s’effectue en repérant 
la teinte dans la région M. Supposons que dans les 
régions planes (loin du contact) lappareil soit régle 
pour avoir le pourpre du premier ordre, par exemple 
(0,565 micron). Si on observe le bleu ciel du 
deuxieme ordre en M (0,664 micron), on aura: 
A — 8 = 0,664 — 0,565 = 0,099 p. On verra plus loin 
que cette mesure peut s’effectuer plus commodément 
4 Vaide de franges. Dans la formule (4) intervien- 
nent également d et l’indice n. Le dédoublement d 
est une caractéristique, mesurée une fois pour toute, 
de V’appareil. I] faut done mesurer indice si celui-ci 
vest pas connu. La mesure peut s’effectuer en meme 
semps et de la méme maniére que celle de A — 8. On 
sntroduit dans le liquide un prisme de verre d’indice V 
sonnu voisin de n et d’angle A (fig. 6). Le chemin 


yptique le long du rayon (1) est, a une constante pres, 
/ 
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e’ est l’6paisseur de prisme traversé par le rayon (1) et 
‘ 8 . . 5 . 
e,, ’épaisseur maximum du prisme. On aura donc 


dA, N de’ ; 
aS ant —n) 4, —(W—n)tg 4; 
da‘ 
ar représente la pente de la surface d’onde aprés 


traversée du prisme. 

Soit A’ la différence de marche mesurée par la 
teinte dans une région quelconque du prisme (la 
teinte est constante puisque la pente est constante) et 
8’ la différence de marche mesurée par la teinte obser- 
vée a coté du prisme. On aura 


DN a ON! 
(VN —n)tg A =- 7 , 


dou 
: A‘ —8' 
(6) nN 


7 d tg An : 


IV. Précision. —Calculons d’abord Jerreur sur 
Pindicen duliquide. Supposons que la différence n — NV 
corresponde a une valeur de A’ — 8’ égale 4 2 2 et 
que l’erreur sur A’ — 8’ ne dépasse pas 2/20. En 
négligeant les erreurs relatives sur d et sur tg A, l’er- 
reur relative sur la différence n — JN est égale 


Anal jn) AT aa 
n—WN nN 
Pour une différence des indices n — NV = 10-1, 
Verreur absolue sur l’indice n sera 
AS [hse NOS 


L’erreur sur ® peut se calculer a partir de (4). L’er- 
reur relative sur la différence A — 8 peut étre prise 
égale 4 1/40 comme précédemment. L’erreur relative 
sur m.— 1 sera 

A(n — 1) An 


a eiy s\n Meee: 
= = eA Ome 
n—i1 0,5 
L’erreur relative sur ® égale a la somme de l’erreur 
relative sur A — 8 et sur n — 1 sera 


n— 1 


AB 1 : 
— = — + 5.10? = 0,03. 
B 40 ci 
En prenant par exemple % = 5°, on aura une 


erreur absolue AB = 9’. 
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V. Description de l’appareillage. — La figure 7 
montre le schéma du dispositif utilisé. Un collimateur 
formé par la fente F et l’objectif O, éciaire la cuve G 
contenant la lame immergée L. La lumiére traverse 
ensuite un polariseur @, le systéme biréfringent S, un 
viseur formé par l’objectif O, et l’oculaire Os, puis un 


| 

| 
Ri hey 

| 

| Se 
<——————_—— > 

| oe 


analysateur .{b. On effectue la mise au point du viseur sur 
le ménisque 4 étudier. L’objectif O, a 40 mm de lon- 
gueur focale et on utilise un oculaire O; de grossisse- 
ment 6. Le systéme biréfringent S est un polariscope 
de Savant : il est constitué par deux lames de quartz 
identiques taillées 4 45° de axe et croisées entre 
elles. L’épaisseur totale du polariscope est de 4,2 mm. 
Puisque les deux lames ont leurs sections principales 
croisées, l’onde ordinaire produite par la premiére 
lame devient extraordinaire dans la seconde et inver- 
sement (fig. 8). 

La figure 8a montre la marche des rayons et les 
figures 8b et 8c représentent les traces des rayons sur 
un écran placé entre les deux lames (fig. 8b) et aprés 
les deux lames (fig. 8c). Le décalage définitif que nous 
avons appelé d est O,E, sur la figure 8c. Il est égal 
a7 X 10-*e; si e est l’épaisseur de chaque lame. Les 
deux ondes émergentes correspondant a une incidence 
normale sont en phase a la sortie du polariscope. II 
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ne produit dans ce cas qu’un décalage latéral. On peut 
produire un déphasage entre les deux ondes dédoublées 
en inclinant trés légerement le polariscope. Ceci est 
obtenu par un dispositif mécanique simple, reliant 
le polariscope a l’objectif Oy. 

Notons qu’il est possible d’employer un prisme d 
Wo taston a la place du polariscope de Savart, 
mais ce prisme devra étre placé au foyer de l’objectif 
O,. 

VI. Emploi des franges a l’infini du polariscope. — 
Nous avons vu au § IlI comment peut s’effectuer une 
mesure par la méthode en teinte plate. On peut opérer 
également par le procédé suivant : intercalons entre 
Poculaire O; et l’analyseur un second polariscope de 
Savart. Ce polariscope n’étant pas éclairé en lumieére 
parallele, projette dans le plan d’observation ses 
franges a l’infini. Dans le cas du polariscope de SAVART 
ce sont pratiquement des droites paralléles et équidis- 
tantes. Ces franges dessinent les lignes de méme diffé- 
rence de marche entre les deux ondes &, et Z,. Les 
déformations de lVonde produites par le ménisque 
provoquent des déplacements de franges. Ces déplace- 
ments mesurés en interfranges donnent les différences 
de marche A et 8 de (4) ou de (6). Ce sont des mesures 
tout a fait analogues aux mesures interférentielles 
ordinaires mais ce sont les dérivées du chemin optique 
qui interviennent et non plus les chemins optiques 
eux-mémes. Notons que l'image observée par le viseur 
étant normalement a l’infini (ceil non accommode) le 
deuxiéme polariscope de Savart ne produit aucun 
dédoublement de limage. 


VII. Observation des gouttes liquides déposées 
sur un plan. — Nous avons précédemment décrit la 
méthode d’observation des gouttes. Si la lame solide 
est transparente, on éclaire l’ensemble de la lame et de 
la goutte de liquide par transparence. L’angle de con- 
tact 6 est donné par (4) (fig. 9). 

A — 8 est déterminé soit par la méthode en teinte 
plate (lumiere blanche), soit par la déformation des 


~ 
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Fig. 9. 


franges a l’infini du second polariscope (lumiére mono- 
chromatique). 

Si la lame solide est opaque et réfléchissante, on 
observera le systeme en lumiere réfléchie. 

La figure 10 montre la photographie d’une goutte 
d’eau posée sur une lame de verre non parfaitement 


Fig. 10. — Diamétre moyen de la goutte d’eau : 


1 mm. 


propre, en utilisant la méthode des franges a l’infini 
du deuxiéme polariscope. Les franges sont dirigées 
dans le sens du dédoublement. 

La figure 14 reproduit schématiquement la disposi- 
tion des franges. Considérons le point A ot la tangente 


\ 
J 
: 


LEI 


Direction Gu 
dédoublement 


produit par le 
‘\“polariscope. 


Franges produites par 
le 2°™ polariscope. 
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au bord du ménisque est perpendiculaire aux franges. 
On a immédiatement l’angle de contact en A. Lors- 
qu’on passe du plan au ménisque, dans la région A, 
la différence de marche varie de 3 2/4 dans le cas de la 
figure. La partie entiére de l’ordre d’interférence est 
déterminée en repérant le déplacement de la frange 
noire en lumiére blanche. On aura donc 


A—s8= = (A = Ou,546). 

L’indice de leau est n = 1,33 et le dédoublement 

d = 14,7 x 10-* mm, d’oti langle de contact 8 : 

6B = 4°50". 
On aurait langle de contact au point B opposé de la 
meéme fagon. 

Pour avoir sur la méme photographie l’angle de 
contact en un point M autre que A ou B, un calcul 
supplémentaire est nécessaire. En effet le dédouble- 
ment qui intervient dans (4) donnant l’angle de con- 
tact 6 n’est plus égal a d. Le dédoublement en M s’effec- 
tue dans la direction MM’, or on veut langle de con- 
tact @ correspondant a la section normale PP’ passant 
par M. Le dédoublement des deux ondes 4, et 2%, 
dans la direction PP’ est donc d cos y ot y est l’angle 
formé par MM’ et la trace PP’ de la section normale. 
L’angle de contact 8 en M est alors donné par 


A —8& 
(n —1) dcos y~ 


(7) = 


La relation (7) montre que les mesures ne sont plus 
possibles aux points Cet D ot la tangente au ménisque 
est parallele a la direction du dédoublement. On peut 
éviter le calcul de y en faisant tourner la goutte avec 
son support sur la platine du systeme d’observation 
de maniére a faire occuper au point M une position 
analogue a celle de A par rapport aux franges ou au 
dédoublement. 


VIII. Angle au contact d’une lame plane immergée. 
— On détermine de la méme maniére la pente de la 
surface liquide et angle de raccordement 6 avec lalame 
(fig. 12). Si la surface liquide se releve notablement au 


Fie. 12. 


contact de la lame solide, observation de l’angle 6 
peut devenir impossible. I] peut en effet y avoir 
réflexion totale dans le liquide, arrét des rayons émer- 
gents par la Jame solide MN, ou bien les rayons émer- 
gents trop déviés ne pénétrent pas dans l’objectif du 


84 M. 
viseur. On utilise alors le tilting-plate de ADams 
Jessop : la lame solide est inclinée de maniére a ce que 


la surface liquide reste plane ou presque plane au 
contact de la lame solide (fig. 13). La mesure est alors 
possible si la lame solide est transparente. 


Fie. 13. 


IX. Profil du ménisque au contact d’une lame 
plane. Forme des franges. — L’équilibre du ménisque, 
équilibre entre les forecs d’adhésion et le poids de 
liquide soulevé s’écrit (fig. 14) : 


a 2 Bee ? 
(8) Ala+a = aca R, to} Y; 


R, et R, sont les rayons de courbure principaux du 
ménisque, y l’ordonnée du point M considéré, A la 


4 


Fig. 14. 


tension superficielle du liquide, D sa densité et g acceé- 
lération de la pesanteur. Dans le cas d’une lame plane 
i =O et 

i) 

Fin “ads 


ou ds est un élément du profil du ménisque. 
On a d’autre part 


OO aS Re has 
iis igs 


d’ot, en remplagant dans |’équation (8), 


0) = COSmO 


dé 
Raa = 
aie Dgy 
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et, en intégrant, 
D gy? =—2 A sin 0+ C%, 


pour 0 = 90°, y = 0, done la constante d’intégration 
est égale a 2 A. 

On aura 

(9) D g.y? = 2A (4 —sin 6). 


C’est Péquation réduite du ménisque. 
Pour connaitre la forme des franges, il faut connaitre 
la différence de marche en fonction de z. Ona 


dx 


= jes = teres 
dy 


=} 

L’équation du ménisque devient 

2A (1 —sin 6) 
Dg 


el, par intégration, 


ee a 
2vV¢ 


eo M2+V 1 + sin? ae 
Be VY 1+ sin 0 


La_ constante oe peut étre déterminée 
si on connait langle de raccordement car, pcur x =0, 
8 = 0 angle de raccordement (fig. 14), d’ot 


ns ie Us aemeRs 
y 5) i eter 
oe V2EV 1 + sin % 
a Cee 


On peut remarquer que l’angle 6 (voir fig. 15) est égal 


Dg 


da cotg 8 = A cos 0 d0 


gy) 


‘4+ sin 6 b. 


419.\/ tein 6 


Tw 
; a : rs 
a 8 et par conséquent l’équation (4) s’écrit ici 


A—8 = (n—1)dtg(4 —0), 


Les équations (10) et (11) permettent de déterminer 
la forme des franges. En donnant a @ une valeur com- 


(11) 


prise entre et == Peet “ Péquation (10) permet de 


calculer Pabscisse x correspondante. La méme valeur 
de 6 remplacée dans l’équation (11) donne A — 8, 
donc le déplacement de la frange exprimé en inter- 
franges. 

Dans le cas du contact eau-plexiglas & 15° on a, 
en enfongant la lame, 


= 73,49 dynes/cm, Sr 76e 
he = 133. NG os 10-3 mm. 
Le tableau suivant donnant les valeurs de A — 8 


en microns et en interfranges en fonction de 2, dis- 
tance a la lame immergée, et de © pente du ménisque. 


. ‘ 
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Q° So A —§ microns | in =P gee u) | 
76° | 0 1 2212 | 
ry OYAT: | OG ees 
ee | Om 1,16 2,132 
| oe 112 ee 
77°36 0,27 1,08 1,972 | 
78°36 | 0,54 | 0,99 1,80 » 
Pe) | 0,67 0,04 1,76 » 
7993, | 0,81 | 0,90 1,64 » 
80° : 0,94 0,85 1,56 
80°30 | 1,02 0,81 1,48 
Snoe | 1,24 | 0,77 1,40 
os 1,41 0,73 | 133 
Soe | 1,44 0,68 1,25 
82 10) 1,68 0,64 | 1G 
830) 1,88 0,60 1,09 
83°30 1,99 OF55 1,01 
pes 2,30 0,51 0,93 
84°30 203 0,47 0,85 
85° 2,74 0,42 0,77 
86° 3,34 0,34 0,62 been 
870 17, 0,28 0,16 bs Sma 


La figure 15 montre la forme calculée et Ja forme mesu- 


Krom 


par la méme relation que précédemment : 


est le diamétre du fil, on a 


5 R page - 
== 96 = des Mics > 
; 2S lt eee 
4 x 7 
Nous avons done 
A-d — 
a Ga 
O O OS 18 sFoyn x A | cos 0 dy aie we ot =-—Dgy 
: a + eae 
IRN 1S). : 5) 


rée des franges et la figure 16 l’aspect des franges 
dans le cas du contact eau-plexiglas. On pourra noter 
les petites variations locales de langle de raccorde- 
ment. 

Notons que les équations (10) et (11) permettent de 
déterminer la tension superficiclle A et angle de rac- 
cordement . On mesure deux valeurs de 0 (équa- 
tion 11) et les deux distances # 4 la lame qui leur cor- 
respondent. On a alors deux équations analogues a 
Péquation (10) qui permettent de calculer A et 6p. 


X. Angles au contact d’un fil cylindrique immergé. 
—L’observation du ménisque d’un fil cylindrique se 
fait avec le méme dispositif que celui qui a servi a 
Pobservation du ménisque d’une lame. Nous allons 
chercher a calculer dans ce nouveau cas particulier la 
courbe A — 8 = f(z) c’est-a-dire la forme des franges. Fie. 17. 


L’équilibre du ménisque se traduit encore par la rela- 
tion (8). Le rayon de courbure principal R, est donné 


mais le rayon de courbure principal Ff, n’est plus 
infini. On peut l’exprimer en fonction de a. Sur la 
figure 17 la partie hachurée représente le fil et lori- 
gine des abscisses est prise au contact avec le fil. Si a 
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et, en tenant compte de la relation 


daz 

tg §@= —— , 
: dy 

on obtient 
dé I 
A |—cos 6— tg 6+ = —D gy. 
dz a 
Th ea 


Dérivons par rapport a a: 
d?6 AON | 
sin 6 —— + cos 6 | ) + | 


da? \da/ / a 


ali 
+ Dg cotg 6 = 0. 


SEIS: 


En intégrant on obtiendra la relation cherchée entre 0 
et x, c’est-a-dire une équation analogue al’équation (10). 
L’équation (12) est une équation différentielle du 
second ordre qui ne peut étre résolue d’une facon 
rigoureuse. Le calcul numérique a été fait par M. Bur- 
TIN que nous remercions vivement. 

En prenant les valeurs numériques suivantes : 


A = 73,49 dynes/cm, 


6, = /6°, — = 0,17 mm. 
0 ) 9 p) 


Le tableau suivant donne les résultats pour quelques 
valeurs de z : 


3) zr + “Oy 
76° 0,100 mm 
77° 0,104 
78° 0,109 
79° 0,114 
80° 0,120 
81° 0,127 
82° 0,136 
83° 0,146 
84° 0,161 
85° 0,181 
86° 0,210 
87° 0,266 
88° : 0,389 
89° 0,684 


A partir de ce tableau, on peut déterminer la forme des 
franges. 

La figure 18 montre la photographie des franges 
obtenues au moyen du dispositif de la figure 7. Ces 
franges ne dessinent pas les lignes d’égale pente de la 
surface liquide, mais les lignes d’égale différence de 
marche A — 8. D’aprés Véquation (7) ces franges 
caractérisent donc les lignes pour lesquelles l’expres- 
sion B(m — 1) d cos y est constante comme dans les 
cas précédemment étudiés. 

Voici comment on peut déterminer la forme des 
franges d’apreés les valeurs numériques du tableau. On 
trace une série de cercles de centre O (fig. 19) et dont 
les rayons sont donnés par les valeurs de X (colonne 
de droite) du tableau précédent. Ces cercles représentent 
les lignes d’égale pente du ménisque entourant le fil 
de diamétre a. Par exemple le point M est sur le cercle 
de rayon 0,27 mm auquel correspond une pente du 


[ Orr, Acta 


tt RT ISS Tr 


Fig. 18. 


ménisque 8 = 3°. Cherchons le point M sur le cercle 
de rayon 0,39 mm et qui soit sur la méme frange que 
M. Soit y,, et 8,, les valeurs de y et 8 qui correspondent 
au point M. Soit y, et 6, les valeurs de ces paramétres 
qui correspondent au point P. L’équation (7) indique 
que les franges sont les lignes pour lesquelles la quan- 
tité B(n — 1) dcos y est constante. Puisque n et d ont 
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des valeurs constantes les franges dessinent les lignes 
le long desquelles le produit 8 cos y est constant. On 
a donc 

Px COS Yue = cm COB Yp> 


dou la valeur y,, fixant la position du point P sur le 
cercle de rayon 0,59 mm, 


M 
C03 == = ——.008 ¥-: 
P 


Le point P se trouvant sur un cercle le long duquel 8, 
est égal a 2°, on aura 


cos 80°30' 


cos ¥,, 


87 


a 


Np sues 


2P 


| queWweajqn 


dou y, = 75°90’. On peut ainsi tracer point par point 
les franges d’égale différence de marche A — 8 qui 
sont tracées sur la figure 19. En projetant un systéme 
de franges données par un deuxieme polariscope on peut 
faire les mesures comme il a été dit au paragraphe VI. 

La figure 20 représente l’aspect des phénoménes 
dans ce dernier cas. Le fil a été incliné 4 45° et appa- 
rait comme une ligne verticale sur la photographie. 
La petite tache centrale indique le point ot le fil 
touche la surface liquide. 


XI. Limite d’emploi de la méthode. — Pour que 
la méthode soit applicable il faut que les rayons lumi- 
neux puissent émerger c’est-a-dire qu’il n’y ait pas 


tate) M. 


réflexion totale dans le liquide. I] faut également que 
les rayons lumineux pénétrent dans!’objectif dusystéme 
dobservation dont l’ouverture doit étre suffisante.La 
premiére condition peut toujours étre facilement réa- 
lisée dans le cas d’une lame plane immergée en ut li- 
sant le procédé du tilting-plate. Toutefois, si la lame 
solide est opaque, ce dernier procédé ne peut pas étre 
employé. La lame étant maintenue alors verticalement 
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le liquide peut se relever notablement au contact de 
la lame et produire la réflexion totale. On peut alors 
opérer comme il a été dit au § 9 en mesurant deux 
valeurs de la différence de marche, c’est-a-dire deux 
valeurs de 6 éloignées de la lame solide. Pour mesurer 
les valeurs de « correspondantes, sans étre géné par la 
réflexion totale, on pourra repérer la lame en l’éclai- 
rant latéralement au moyen d’une lampe auxiliaire. 


Manuscrit regu le 21 février 1958. 


Determination of the intensity distribution resulting from the random 
illumination of a plane finite surface 


Lt. Marx J. BERAN, 


Antenna Laboratory, Electronics Research Directorate 
Air Force Cambridge Research Center 
Air Research and Development Command Bedford, Massachusetts 


Summary. — In this paper the intensity distribution resulting from random illumination of a plane finite surface is found by solving 
the equations governing the propagation of the two point correlation function in free space. Using an ensemble average the dis- 
tribution at some distance from the radiating surface is found to be the LAMBERTIAN distrib:tion multiplied by the FouURTER 
Transform of the correlation function defined over the illuminated surface. The ensemble average is then shown to be equivalent 


to either space or time averages in a variety of cases. 


SomMAIRE. — Dans cet article, on calcule la distribution de Vintensilé provenant d’une illumination quelconque d’une surface plane 
finie, a partir des équations qui régissent la propagation de la fonction de corrélation de deux points de Vespace libre. Utilisant une 
moyenne d’ensemble de la distribution a quelque distance de la surface éclairée, on trouve que la distribution est celle de LAMBERT 
multipliée par la transformée de Fourier de la fonction de corrélation définie sur la surface éclairée. On montre ensuite que dans 
toutes sortes de cas la moyenne d’ensemble est équivalente aux moyennes dans l’espace ou dans le temps. 


ZUSAMMENFASSUNG. — In dieser Arbeit wird die Intensitdtsverteilung berechnet, die sich aus der ungleichmdssigen Beleuchtung 
einer ebenen begrenzten Platte ergibt. Das geschieht durch die Losung der Gleichungen, die die Ausbreitung der Korrelations- 
funktion fiir zwei Punkte im freien Raum beherrschen. Indem man ein Gesamtmittel der Verteilung in einiger Entfernung von 
der strahlenden Fldche benutzt, erhdlt man fiir die Verteilung das LAMBERTSCHE Gesetz multipliziert mit der Fourtertransfor- 
mation der Korrelationsfunktion iiber die beleuchtete Fldche. Das Gesamtmittel ist dann in einer Reihe von Fdllen gleichwertig 


mit dem rdumlichen oder mit dem zeitlichen Mittel. 


1. — When a plane finite surface is illuminated in a 
random manner the assumption is often made that 
the far field intensity distribution is Lambertian. 
This is correct if the variations in field on the illumi- 
nated surface occur in distances of the order of a wave 
length or less. However when the field variations, 
though random, are not this rapid, the intensity 
distribution is not Lambertian. In this paper it is shown 
that in general the intensity distribution is the Lam- 
bertian distribution multiplied by the Fourrer trans- 
form of a correlation function defined on the illumi- 
nated surface. 

We shall begin by giving a general definition to 
random radiation. The definition may not appear to 
the reader to be appropriate for the problem discus- 
sed above but it will be shown in section II to be 
applicable for a variety of illuminating mechanisms. 

The definition used will be an ensemble definition 
and it will be for a scalar variable, the latter simplifi- 
cation in optics being supported by the work of 
Green and Wo r [1]. We consider in free space a 
finite surface S enclosing a volume V and suppose a 
very large number of experiments will be performed 


by impressing on S in the nt experiment a scalar field 
of the form : 


NP E186) aL A(X vovneode 


—> 

(where X, represents a coordinate with three posi- 
tional rectangular components x, y, and z and the 
subscript s indicates coordinates on S) and observing 


in each experiment the field " &(X, t) at all points 
external to S. 

It will be assumed that in any experiment our a 
priori knowledge of the field on S is confined to corre- 
lation functions of the form : 


N 


ee Sts ; 1 —> yb wi SES 
(2) J (Xai XsJ= ie Ne » "E,( X30) B*(Xo,0)e 
n=0 


‘EA Xe «) being the complex conjugate of “F(X on 


; —> => 
* It is noted that J..(X.;, Xg,) is in effect the definition 


of the unnormalized complex degree of coherence between two 
points j and k on a surface S. 


OP lrhieA- AGEA 


Volume 5 — 1958 


ADMINISTRATION : OPTICA ACTA, 3 et 5, BOULEVARD PASTEUR, 
PARIS" 15¢ 


BOARD OF CONTROL — COMITE DIRECTEUR — KURATORIUM 


J. M. Orero y Navascues (Madrid), Président. 

R. W. Dircupurn (Reading, G. B.), Secrétaire. 

A. ArnutF (Paris), A. Brot (Gand), J. Casannes (Paris), P. FLeury (Paris), 
F. Gasier (Vienne), G. Hansen (Oberkochen), E. Incersram (Stockholm), 
H. Korre (Braunschweig), W. S. Stites (London), G. ToRALpo pi FRANCIA 
(Florence), A. C.S. van HEEt (Delft). 


EDITORIAL BOARD — REDACTION — SCHRIFTLEITUNG 


English editor : C. G. Wynne, 13, Elwill Way, Beckenham, Kent. 
Rédacteur francais : A. MARECHAL, 3, boulevard Pasteur, Paris, 15°. 
Deutscher Herausgeber : G. Franke, Laufdorfer Weg 2, Weitzlar. 


PUBLISHER — EDITEUR — VERLAG 


Société de la Revue d'Optique, 3, boulevard Pasteur, Paris, 15°. 


Abonnement au tome 5: 3500 F ou $ 10 
(environ 192 p.) valeur sur Paris 
a adresser a la Revue d’Optique, 3, boulevard Pasteur, Paris 15°. 
Conditions spéciales pour les membres des Comités nationaux d’optique 
(s’adresser a ces Comités) 


CONTENTS — TABLE DES MATIERES — INHALT 


M. J. Beran: Determination of the intensity distribution resulting from the random illumination of a plane 
finite surface 


OER tee ir 5 sot PE ane Se Tre ONE aa Rae RA eS Aa Or gM a ne SE And tt ee ee 64, 100, 


D. Canars-Frau et Mme M. Rousseau : Influence de l’éclairage partiellement cohérent sur la formation des 
maces de quelques objets étentlus opaques’..-1e% exh a: 4. 6.50 ee ee een eee 


Bibliographie 


R. S. Crisp : voir P. FisuHer. 


P. Fisner, R. S. Crisp and S. E. Witiiams : A photon-counting grazing incidence spectrometer for the 50- 
AN OOOBAR AT Semek OD  ceohila ale oh, 0, ROKR ON, | tae ee, Cer tc cee: tee 


M. Francon et Mle Y. Gannon : Etude des angles de raccordement des surfaces liquides avec les solides par 
losmnterierences. ery lumibre, polarisee .... ass acne sesame ae ae ene eee 


Y. Ganpon (Me) : voir M. Francon. 

R. G. GiovaNe ii : The influence of scattering within photographic emulsions on resolving power.......... 
RCT OL CONS Man tee Peer ae chk RY INET, «ile aaah Uys gt sh, AS aS gg 63, 
F. D. Kaun : On photon coincidences and HAnBuRY Brown’s interferometer ..............00 00.00 eee eues 
peCcelcELLY = lhe interferometry measurement.of metal oxide films... ic cnuslisu wis awisbsne os pee te Se 
Poti aLInvoot-Qualhty evaluations ol opticalisystems ... 24.5 oi Uiee Ae ec eee ee 
A. MarEcHAL : La diffusion résiduelle des surfaces polies et des réseaux... . 0.00.0... cece cee ee eee aes 
H. Marx: Linearisierung der Durchrechnungsformeln fiir windschiefe Strahlen (I])...................... 
M. RoussgEau (M™e): voir D. CaANAts-FRAU. 


O. StierstapT und W. Zwiks.ter : Einige Reobachtungen uber die Wahrnehmbarkeit kurzdauernder, farbiger 
Eichterschemungen und, ihre: Deutung..:......v5..cceh 2 sete cos Sis dg muee & noegl e ciee : SOU ern ce 


ieelsuotucHt: voir Y. UKITA- 

NMeWxira and J. Lsusrucui : Interferometric study of microscope objectives ...........4.24.s.s meet een 
A.C. S. van Hert and A. WALTHER : Some methods for measurements on thin films...................... 
A. WautHeErR: voir A. C. S. vAN HEEL. 


FE. Wasser : Investigation into the theory of prediction of the appearance of colours and its bearing on the 
PHCOTY (Ole COLOUL “VISION y. 5. wftte o sfase orm. 9 Pics, aoe nar nOLk bu iotAte Se sales ces ani ay enn es see eee 


S. E. WiLuiaMs: voir P. FISHER. 


W. ZwiESLER: voir O. STIERSTADT. 


Pages 


88 
109 


ANALYTIC SUBJECT INDEX — TABLE ANALYTIQUE DES SUJETS — SACHVERZEICHNIS 


Beugung (s. u. Diffraction). 


Bibliographie (voir Book Review). 


BOOKSTCVIC WHE ee te Ree ase eae 64, 100, 
Biicher (s. 0. Book Review). 


Colloquium 


Controle (voir Testing). 


Couches minces (voir Thin layers). 


Diffraction: Influence de l’éclairage partielle- 
ment cohérent sur la formation des images. 
Intensity distribution resulting from the ran- 
dom illumination of a plate finite surface. . 


Diffusion: The influence of scattering within 

photographieremulsions =, 45.3. seca 
Diffusion résiduelle des surfaces polies et des 
reseaux 


Diinne Schichten (s. u. Thin layers). 


Geometrical optics : Linearisierung der Durchre- 
chnungsformeln fiir windschiefe Strahlen. . 


Interferometry (see also Testing and Thin layers). 
On photon coincidences and Hanspury 
Brown’s interferometer 


Pages 


109 


93 


Paves 


Information theory : Quality evaluations of opti- 
cal systems 


Instruments: A photon counting spectrometer 
for'50-4 000 Avranges acne eee 

Interferenz- Verfahren (s. u. Interferometry). 

Oplique géométrique (voir Geometrical optics). 

Priifverfahren (s. u. Testing). 

Sehen (s. u. Vision). 


Spectrometry : (see also Instruments). 


Testing (see also Thin layers). 
Interferometric study of microscope objec- 
tIVOS iseutte Le ME e Ae ak ad oe aoe 
Etude des angles de raccordement des sur- 
faces liquides avec les solides par interfé- 
rences 


Thin layers (see also Testing). 
Some methods for measurements on thin 
fils} x 4 ES 
Interferometric 
oxide films 


measurements of metal 


Vision : Wahrnehmbarkeit kurz dauernder, farbi- 
ger Lichterschemungen "27.0.... 0 eee 
Prediction of the appearance of colours... . 
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— => . . . 
J (Xj, X,,) is defined similarly for points external 
10) Sh 
It is in the sense of the above defined correlation 
functions that the field is random. The only meaning- 
ful manner in which we shall discuss the field 1 in the 


above context is to determine how De i) is 
determined ou be Oe Xe ,): It is noted that the 
intensity 1 s(X;) resulting from the carrier frequency 
®, 1s simply J Oe x, 


In free space the scalar field ’ 'E(X, ) is governed by 
the scalar wave equation : 


(3) Tae HOCeayn eK SEX @) 80 


Oe 


To solve equation (3)"F ®) must be given over 


Ss 


the boundary and OBEX. w) must satisfy the radiation 
condition at infinity. 


To find Jes, O) from a knowledge of / se iG i) 
use must be made of equation (3). This may be ‘lone 
as follows : : Multiply equation (3) by the constant quan- 
tity "E,*( X,, 0 ). This quantity may be brought inside 


the differentiation since it is a constant. Hence 


(4) VY [BOS yee. «(x 5p ae 
Searle He X, «) ave *(X. po] =o. 
Next average over n experiments asn—>o. The 


average of the preduct "E(X, S): ee. (X50) | is the quan- 
tity defined above as J C(x se ;): The av eraging over 
all experiments may be pertormed within the differen- 
tiation since sufficient regularity of the functions is 
assumed to allow for the interchanging of infinite 
summation and differentiation. Thus we have 


(9) ie J Xs X.,;) si K* J A(X; X,;) ear (7) 


To solve this equation it is only necessary to know 
J RX x. j), that is the correlation between the point 
xX, and all other points on S. The solution then repre- 


ae the correlation between the point X; and all 


other points in space. 

What is really desired however is the correlation 
between any two points in space, Let us call these 
points xX; and en Setting X= X; in the solution of 
equation (¢ (5) aves the correlation between the points 
cK and X.. If now the equation is solved with X,; as 
a variable (though fixed in each solution) the corre- 


‘ation between x; and all Pos on the surface may be 
“etermined since J AG: Be J=/d U(X, ie j): 
() A similar result for a time averaged correlation function 


ras obtained by E. Wowr (Pro*. Koy. Soc. A, 280, 1955, 
p. 246-265). 
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To find TK xX) the equation 


(6) WJ(%,%) + K(X, X) =0 


is then solved for J eee x) using as a boundary 
condition J ES. X)). 


The above croeedGre shows that the correlation 
function of two points external to S may be determi- 
ned from a knowledge of the correlation of all pairs 
of points on S. Since NAGS) = J (Xj, X;) this means 
that the intensity distribution in space may be found 
if the correlation between all pairs of points on S 
is known. It must be remembered though that 


di Oe, De ‘) is complex and that in general both its 
real and tice ners: parts must be known over S. 
Before interpreting J,,(X,, X;) in terms of the 


original statement of the problem let us find TCO) 


from a knowledge of ener x) when S is a plane 
finite surface. 

We first set up a coordinate system with the xy 
coordinate axes in the plane containing the radiating 
section and the z axis perpendicular to the plane. 
The origin will be located at some central point within 
the radiating section. The area of the radiating sec- 
tion will be denoted as A with a characteristic dimen- 
sion a. 

The boundary condition will be that the field is 0 
on the plane z = 0 except over the area A. The justi- 
fication for the discontinuity is that 1/K <a. 

To characterize the correlation function of the field 
on A, we define |/| as a characteristic length over 
which the correlation extends. When J, (oe $¢ ;) 18 
integrable / may be defined as 

MEP rad , Te OC 5 € dX 
=gelbe | sk.“ ” sk ° 
It is noted that / is real when the statistics in some 
surface or region are homogeneous and isotropic. Ho- 
wever whether or not J eG ~ NE. 4) is integrable we 
assume it approaches 0 as x, ~~ Ser > ©. 

The problem will be solved under the assumptions 
that a/l > = 1. Ka>|, | Kz > 1, zjas>1, Kal/z<1 and 
Fo hp weX ee a ge BD 
The ies booadition 1S eee as meaning that 
the correlation function depends only upon the dis- 
tance between two points on S and not upon their 


absolute positions. The imaginary part of J mee 6 i) 


is therefore O on S. 
To begin we solve equation (5). Using GREEN’s theo- 


rem, we find that an integral solution for EG X,) 
isc: 
Hees MOG ae 


a Ke = lkh 
Ae 5 ee = X, |) Re da, dy,’ 


90 
where R = [ (a —a,")? + (y—y,/)? + 22 ]¥? and K has 
been assumed to be > 1/R. 

Next expand R ina binomial series, transforming to 
the coordinates 2} = 2, — i, Ys = Yor — Ysi and 
x! =2—4,,y' =y¥—Y,j- Pines KlP/z<1 and x} and 
y, are of order 1, R may be approximated in the ‘term 
exp (i KR) a 


Bale en ey; 

: — ala’ —yly 
aha 

Ela oa Me 5 


In the term 1/R? only the first term is necessary, so 
that the approximate expression for J,,(X, X;) is 


eK eiK(a'2-+y/2+22)8 


(9) Oe Xa) = ; Z 
Ay Aa) = OR Gey? be) 
i K (ego! + yy") 
x ff 3 (|X?) e @t+yt ait da’ dy, 


If now X is considered a particular point ea 
eds : —> we cee 
Xs; is replaced by the general point X,, J,(X,, X)) 
i Kz; e—ik| (xj—a,)? +(yj—ys)? +23 
De ae 

m [(%,—2))*+ (yj —y,)?-+ 25] 

ik [xt (x; —2,) se ysl; —ys)] 


*(|X"]) e [i —)? + (yj ys)? + 25] da’, dy! 


(LOVE (a Kas 


se festa 


since Tee ease. (Xs dite 
This expression for TA X;) serves as a boundary 
condition for the equation determining the correlation 
. 7 => _ 
function J..(X, X;), the latter function being also 


governed by equation (5). Using a GrreEn’s function 
(assuming KR > 1) 


DIE (X exe 
Ae ek > > eiKR 
where R =[(x —2,)? + (y—y,)? +27} 


X;), X = X;. Using 


Since we desire here Tees 
(X,, xX) simplifies to 


equation (10) the expression for J, 


K2 Ze A 
An? (a? + y3-+22)? 


i) we LK (e' ZX. =e y! “Ves 
{pi Jo*(| XG) exp = ett d 
where the condition Ka l/z<1 has been used to approx- 


é iKu X 
imate such terms as exp 


s 
eee 7, 
(+ y+) , 


(12) J (Xj, X,) = 


ay dy! 


s 
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When K/ < 1 equation (12) reduces to the lamber- 


tian Aetebotea : 
IAP 
An? (2% 


It is remembered that I(X;) i is aes to Je 


(ye Jp ea aie 


jp 


tl. —It is the purpose of this section to show that 
ti ROG as obtained from equation (13) above can be 
used to determine the intensity distribution (subject to 
the stated assumptions) for the following three types 
of radiation. 

a) Monochromatic light passing through a rough 
surface (such as ground glass). 

b) Unmodulated polychromatic light passing through 
a rough surface. 

c) Randomly modulated polychromatic light pas- 
sing through a rough or smooth surface. 


Il a. —If the monochromatic light (@ = wp») pas- 
ses through a medium such as ground glass it emerges 


in an exit plane with a field distribution that may be ~ 


described as random. Denote the scalar component 
discriptive of the field in this plane as 


=e OO 


HOES) en: 


The quantity 

pee ew ie wee te 

ok, B+ Q=a fa 
Ak 


Al = SF et —> 
af LX 4 @) EA(X, -L oes ©.) dx, 
A 


af 


may be considered the correlation function which — 


describes the randomness, where the integration is 
taken only over the illuminated area of the plane. 


pe ee ee ; 
Since it is assumed that o(X,, X, + &,) falls to zero — 


in a distance J, small compared to a characteristic 
dimension, a, of the illuminated area, edge effects 
are unimportant. It is further assumed that within 


the illuminated area, o(X,, Le Lah By depends only 
al 

upon |&,| not upon the absolute position of X or 

direction of & 

geneous and isotropic. 

The fact that 1 is very small compared to a also 
allows us to equate o(X,, ie ++ re )and J (Xe Le + es ) 
since due to the statistical Romosonmne over A many 
experiments are essentially being performed simul- 
taneously over A. Averaging over all pairs of points 
in A yields the same result as the averaging over a 
particular pair of points during many experiments ; 


: 1. e. the statistics over A are homo- — 


f 


this equivalence being a rough statement of an ergo- 


dic theorem. 

It follows from this that in the terminology used, 
the intensity distribution over A would be constant. 
Actually if it were possible to measure the point 
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intensity at points a distance 1 apart, this intensity 
would vary considerably from point to point. Howe- 
ver when we speak of the intensity here, we mean 
the point intensity averaged over areas Jarge com- 
pared to /*. It is only the variations that are very 
very slow compared to the optical wave lengths that 
are usually considered in practice and it is this space 
averaged intensity that is treated here. 

The problem is not as simply resolved when we 
wish to find a similar averaged intensity at points 
external to A. We can then no longer average over a 
homogeneous surface such as A to obtain an average 
space intensity. Away from A the form of the irregu- 
larities change continuously and in fact the very far 
field may assume a form independent of 1 on A. 

Fortunately however we would expect to find 
regions of what may be called local homogeneity since 
the irregularities should have the same statistical 
form over regions large compared to the dimensions of 
any one irregularity. To be sure, the statistical pro- 
perties change with distance but slowly enough to 
allow their form to be determined in any region by 
some method of averaging. Further these statistical 
parameters, in which we have in mind of course the 
space averaged intensity, should depend only upon 

> —> 
o(X,, X, + &) over A and not upon the detailed 
structure of the field at A. 

Now it is virtually impossible to actually prove 
theoretically by non statistical arguments that the 
field is actually as described above and the moment 
statistical arguments are used we are faced with the 
problem of showing that the statiscal techniques used 
correspond to the physical problem. The position taken 
here will be to assume this description is sound, show 
from this assumption that an ensemble average is 
equivalent to a local space average and then briefly 
cite some evidence obtained by using the ensemble 
approach which bears out the original assumption. 

To equate the ensemble and local space averages we 
note that given a particular random field over A the 
detailed form of the local irregularities external to A 
are determined. However by assumption, in some region 
large compared to the size of a fluctuation the statis- 
tical properties are homogeneous and space averages 
over this region may be used to determine them. By 
hypothesis, these statistical properties depend only 


upon 9(X,, X,+ &,), itself a space averaged statistical 
property that has been equated to an ensemble 
average, and the external local space averaged quan- 


‘ities may be equated to ensemble averages in the 
> 


‘same manner as o(X,, se + &.) has been equated to 
> —> —> 
SF (Xop x, v Es): 

1..(X;) is therefore physically interpreted as ave- 
rege intensity, the average being a space average over 
a local region in which the point intensity fluctuates 
ripidly. These fluctuations will be fast compared to 


ae a . . . % 
the manner in which [,,(X;) will vary once it 1s so 
- determined. 
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The ensemble approach as explained in section I 

shows that the average intensity external to A depends 

> > —> 

only upon the function J.,(X,, X, + &) which is 
consistant with the latter part of the assumption made. 
To determine whether there are regions of local 
homogeneity LEG x, + a must be found and the 
characteristic dimension, l, of its extent, determined. 
This dimension must be very small compared to 
distances in which the expression in equation 12 va- 
ries appreciably. In a work to be published shortly 
the author has shown that in the case when / K 1, in 
the far field this is true. It has also been shown to be 
true for all ratios of 1 K for radiation from an infinite 
plane and for a case of radiation from a sphere. 

On this evidence it is concluded that the assumption 
of local homogeneity external to A is consistant and 
that subject to experimental verification equation 
12 may be used to find the intensity distribution 
resulting from monochromatic light passing through a 
rough surface. 


II 6. — By unmodulated polychromatic light is 
meant here light with a frequency spectrum contained 
in the visible band (@, << ,). To analyse this 
case each infinitesmal frequency band between w and 
o + de is first treated separately and T(X) obtai- 
ned as in the previous section. 

1(X;), the intensity obtained from the entire fre- 


quency band is here defined as 
; 4 (FF > > 
im —— e( Xs) e*( Xt) dix 
T>o ah 0 
This is the usual definition of intensity at visible fre- 
quencies and was obtained in the monochromatic case 
—> > 
without explicit integration since e(X, t) e*(X, Ace 
a i is —j Rw ‘ 
E(X,0) &°! -£%(X,0) € @f = E(X, w) E*(X, 0). 


Substituting for 


=> —>- POs —> . 
e(X,t), e(X,#) = | E(X, 0) & do 
e/ 
yields : 
| ft! 
(14) T(X.) =) lim +f. dia 
T> Les 

PA Ws > iot op) => —io't 

x | E(X, 0) dda | “B*(X, «') edo’. 
J @, 4 


Assuming that the limits and integrals can be 
interchanged : 
—> Oe > 
(15) I(X;) =f dw E(X,) Xx 


e’ 1 
Fr We / 4 tel fe ; \ 
be E*(X, @') dw’ | lim +f eltfa—o’) gy iz 
/ oO; Nien gla 


The latter time integration is a delta function 


8(o — o’) and hence : 
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We —> > 7 We > 
i do E(X;, a) E*(X;, ) =| I (X;) de 
Or 


/ @, 


ee . . ¢ 
where J,,(X;) is determined from equation (12). 
If for some reason the rough surface caused each 
frequency component to change its phase and ampli- 


tude in a different manner JX(| x) would simply be 
a function of frequency and the same formulas would 
apply. 

II c. — This last case of randomly modulated 
polychromatic light is meant to treat the case of vi- 
sible light which is slowly modulated in time by some 
random extended source such as an incandescent bu b. 
Itis supposed that the time modulation is different over 
different parts of the source and that the characte- 
ristic modulation time constant, Jj, is very large 
compared to 1/a@) = Jj) Where w) is the average 
radial frequency of visible light. 

Consider first this light passing through a smooth 
aperture. In this case at any instant the field on S 
looks very much like the field which has passed through 
a rough surface and the equating of space and en- 
semble averaging may again be used to find J (Xj, X,,). 
However it is instructive to show how time and 
ensemble averages may be equated to yield the same 
results without the necessity of assuming local homo- 
geneity external to A. It should be mentioned however 
that the two comparisons are only possible in the 
same experiment when the statistics over S are homo- 
geneous. 


5 > 1 each frequency component of the car- 
Jeo 

rier wave spectrum passes through many cycles before 
the modulation is effective. Let us choose a time 
P such that J. <2<J. Inthe time interval from ¢ to 
t+ T the modulation has negligible effect and we may 
suppose that instead of a continuous process we have 
a discontinuous one in which at times ¢t, t+ 7..., 
t+nf...,t+ NT an experiment is performed by im- 
posing on S a field 

Oe n y iwt 
[ E(X, @) e dw where a, Ca <a, 
Oy 


« 
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contains the visible light spectrum. Since the modula- 


—> 
tion is assumed random, "£,(X, ») has some random 
distribution over the smooth aperture S at time 
t+ nf. 


To find the correlation between two points ue and 


GC for the frequency component , we simply com- | 
pute the sum 


N 
>) "£(X,,@, nt) "E*(X,, , nf) . 


The only difference between this formula and the 
ensemble formulation given in section I is that here— 
"Ff and ™*1E are closely correlated whereas they 
were implicitly assumed to be independent in section I. 
However as NZ is large compared to J the two for- 
mulations yield the same results since the time corre- 
lation is then unimportant in the ensemble average. 

Thus J Lae we .) is determined by appropriately 
interpreting the ate process as an ensemble of expe-_ 


riments. Knowing J, she, X;) we ey then use the 


reasoning in section II b to find Jie x 
When the aperture is a rough puntaes and the car- — 
rier frequencies are slowly, modulated we compute 


separately J, hoe jee the correlation function for — 


sj? 


the rough surface alone and J OS D6 im the corre- 


sj? 
lation function for the modulation on a smooth sur- 
face. These correlation functions should be indepen- 
dent of each other and hence : 


J Oe ie = J nee! ee + J ol, i) Xe im Where 
HERD Gs ber a Lees ie iM a 
be found as in section I b. 
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PHOTON COINCIDENCES AND HANBURY BROWN’S INTERFEROMETER 93 


On photon coincidences and Hanbury Brown’s interferometer 


F. D. Kaun 
Astronomy Department, The University, Manchester 13. 


SUMMARY. — Some statistical results are established 


which govern the rate of detection of photons reaching a small receiver 


and coming from a source which is too small for the receiver to resolve. The results are used to find the probability of the 
coincident arrival of a pair of photons at two small receivers trained on the same source. In essence, HANBURY BROWN’S 
interferometer is a device which measures this probability. In § III the statistical results are used to estimate what signal : 


noise ratio may be expected when the interferometer is used. 


SomMaAtIrE. — On établit quelques résultats statistiques qui régissent la capacité de dtection de photons arrivant sur un petit récepteur 
et issus d’une source trop petite pour étre résolue par le récepteur. Les résultats sont utilisés pour trouver la probabilité d’arrivée 
simultanée, pour un couple de photons, sur deux récepteurs dirigés vers la méme source. L’interférométre de HANBURY BROWN 
est essentiellement un syst¢me qui mesure cette probabilité. Dans le paragraphe III, les résultats statistiques sont utilisés pour 
estimer le rapport signal/bruit qui peut 4tre prévu lorsqu’on utilise Vinterférométre. 


ZUSAMMENFASSUNG. — Es werden einige statistische Ergebnisse mitgeteilt, die die Zahl der nachgewiesenen Photonen bestimmen, 
welche einen kleinen Empfdnger erreichen und von einer Lichtquelle ausgehen, die zu klein ist, um von dem Empfédnger aufgelést 
zu werden. Die Ergebnisse werden dazu benutzt, die Wahrscheinlichkeit dafiir aufzusuchen, dass ein Photonenpaar gleichzeitig 
auf zwei kleine Empfanger fdallt, die auf dieselbe Lichtquelle gerichtet sind. Das HaNpury-Brownsche Interferometer ist im 
wesentlichen eine Anordnung, die diese Wahrscheinlichkeit misst. Im dritten Abschnitt werden die statisiischen Ergebnisse fiir 
die Abschdtzung des Verhdltnisses zwischen Signal und Rauschen herangezogen, welches man beim Gebrauch des Interferometers 


zu erwarten hat. 


I. Introduction. — The new intensity interferome- 
ter invented by Hanpury Brown and T’wiss has 
already been described and discussed elsewhere by a 
mumober of authors ({1, 4, 6, 7, 8, 9, 11, 13, 14}. 
Briefly, it depends on the detection of the coincident 
arrival of quanta at two receivers P and Q, each 
of which is too small to resolve, by itself, the source 
of light under observation. It may be shown that 
the probability of the coincident arrival of a pair of 
quanta is proportional to (1 + ¢), where ¢ is the visi- 
bility of the fringes which would be produced if light 
from P were allowed to interfere with light from Q, 
the path difference being zero. Measurement of this 
probability therefore gives information about .», 
which in turn gives information about the shape and 
the size of and the intensity variation across the 
source of the light. 

Itisimportant to know beforehand what signal : noise 
ratio to expect in such an experiment. The ratio will 
tend to be unfavourable when the intensity of the 
source and the quantum efficiency of the receiver are 
low, and when the time-constant of the measuring 
apparatus is long in comparison with the coherence 
time of the radiation. A rough estimate of some of 
the relevant factors may be made by means of the 
uncertainty principle. As Hanspury Brown and 
Twiss [ 6 ] show, in their § 4, the object of the experi- 
ment is essentially to measure systematic differences 
in phase ®,, = ®, — ®, between observations of 
the classical electric field at P and at Q. Now photons 
are taken from the wave-trains at P and Q, and the 
evatistics of wave-trains are such that even when some 
rhotons have been taken from a particular train, the 
expected number left in the train is still of the order of 
(.©— 1)”. (Here « = hv)/kO, vo is the mean fre- 
euency of the radiation, © the temperature of the 


source). The uncertainty in the number associated 
with the receiver at P or at Q is An, which is of the 


order of n= xQS/2(e® — 1). 


Here 
* = quantum efficiency of the receiver, 
S = equivalent area of the receiver, 


Q = solid angle subtended by the source, 
and 


Ao = mean wavelength of the incident radiation. 
By the uncertainty principle, however [ 10 ] 
An A® = 1 
so that 
rwr(e& — 
Nye 
xOS 
and 
aee—1t) 
MO ase So=<s. 


In the case of most receivers xQS/2¢ is small, of 
the order of 10-4 or 10-5, say. For waves from a radio 
source, ¢ = hv,/k® is of the order of 10-8, or less, so 
that n ~~ 108 or 104%. But in the case of light of, say, 
5 000 A wavelength, coming from a star of 12 000 °K. 
surface temperature, we have « = 2.4, and then 
n= 10-5 or 10-*. Thus 

A®,, ~ 10° or 10-*, for a typical radio source, 


and © 
~z 10° or 10%, for typical visual stars. 


A single measurement on a wave-train from a radio 
star gives the phase difference to within a few minutes 
of arc. In the case of a visual star, however, the uncer- 
tainty is enormous, and even a precision in the phase 
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measurement to, say, 0.1 radian is only achieved 
after sampling 107? or 1014 wave-trains. 

This is the distinction between observations at 
radio and at visual frequencies. One would, in any 
case, use classical theory to discuss the details of the 
experiment at low frequencies and would see no 
reason why it should not be possible. Ina quantum 
picture it is the remaining uncertainty of the number 
of photons in each wave-train which makes the expe- 
riment work. The uncertainty is large in the radio 
experiment and small in the light experiment, with 
converse effects on the uncertainty in phase. 


IJ. Auxiliary Results. — This section deals with 
some preliminary aspects of the problem. We consider 
only sources which have the following properties : 

(i) They emit randomly polarized light ; 

(ii) They are stationary in the sense that the expected 
value of their measured intensity at any instant ¢ 
does not depend on the particular value of t, even if a 
colour filter with any given spectral response is inter- 
posed between the source and the detector. 

We shall be concerned only with sources which 
subtend a small angle, and may therefore divide up 
the radiation field into two components polarized in 
two perpendicular directions normal to the mean 
direction of travel of the light. 

We make a Fourrer decomposition of one of these 
components 


(1) E(t) =f G(v) e— 27 dy . 


If a colour filter is interposed, the component at 
frequency v will, in general, be reduced in amplitude 
by a factor |a(v)| and retarded in phase by an amount 
e(v). The field transmitted by the filter is 


(2) EG) = [ a(v) §(v) e 27 dy , 


where a(v) = |a(v)| el) The expected intensity at 
time ¢ will be proportional to the expectation value of 
the square modulus E(t) 


(3) ([E’@)|?, = 


[ a(v) a*(v') ( &(v) G*(v') ) e2rilv—v")t dy dv’. 
0 


e 7 0 


The value of the double integral must be indepen- 
dent of ¢ whatever function a(v) happens to be cho- 
sen. This is possible only if 


(4) (6(v) &*(v')) = (u(v)) 8(v—v‘) 


where 8(z) is the Drrac delta function with the pro- 
perties that 3(z) = O unless = 0, and 


fe 


| 3(a) da =a4_. 


Equation (4) shows that there are no phase relations 
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between components with different frequencies. We 
now consider some statistical results which follow 
from this property. 

According to equations (1) and (4) the complex 
quantity E(t) is the sum of a large number of other 
complex quantities between which there is no phase 
relation, and the sum of whose mean square moduli is 


J ( 8(v) 8*(v)) dv = (U), say. 


The statistical distribution of the values of 
U(t) = E(t) E*(t) 
is therefore the same as in the random walk problem 


| which is discussed in detail by, for instance, CHAN- 
DRASEKHAR [2] { and we find that 


2 
(U; 


is the probability that |#| shall lie between / and 
L+ dl, or that 


(5) m(l) dl = e PAU) J al 


(6) P(U) dU 1 --UKY) ay 
(U) 
is the probability that | £ |? shall lie between U and 
U + dU. It follows that 


(7) (0) = [0 pO) av = 2( UF. 


It is more useful, for later work, to express this 
result in terms of wavetrains. Suppose the radiation 
field has uniform intensity in the range of frequencies 
Av 

2 


outside. This is a further restriction, but it can be 
removed later without any difficulty. Then &(v) 
may be expressed as a FouRIER series 


Ay : f 
from vy — to vo + ae and has zero intensity 


(8) &(V -- 2) — », E : e2 ire/Av. 


r=— @ 


Now the radiation field will be adequately represen- 
ted if 


( &(v) 6¥(v’) } = ( u(y) ) &v — v’) = up 8(v — v’/), 
since ( u(v) ) has a constant value in the range of 


frequencies concerned. In terms of €, this condi- 
tion requires that 


(9) y » E.€ a e2Ti(re—r'o")/Ay \ =Uy 8(9 — 9’) 
o—ea / 

which holds if (€ ,€ *,) = 0 when rr’ 
Lene oleatd a Aye == rtm 


This result follows because 


and 


* / 
y eee )/Av Ay =: (0 —o') . 
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We have then that 
Be yee 
0 


; Ay/2 
_ sane Give =e ?) ee mivt dy 
—Ay/2 


3 a E_ e27iv(t—r/Ay) do 
—Ay/2 
sin ™ Av(t — r/Ay) 

x Av(t —r/Av) _ 


(414) = e 2 Tivot 


r=>—o 


@ 
— e—27ivot DS €. 


r=—o 


= e-27i vot ‘y € , sine } m Av(t — r/Av) 


r=—_ 0 


This formula expresses the radiation field in terms 
of the well known sinc functions of SHANNON [12]. 
It is seen that the r since function vanishes when 
t = s/Av, where s is an integer other than 7, that it 
equals unity when t = r Ay, and that 


sinc m Av(t —r/Av) — 0 


for large values of | 7 Av(t —r/Av) |. Hence E(t) = € y, 
when ¢t = r/Ay, and this is the instant when the r‘? 
component makes its largest contribution ; the length 
of time during which this component makes an effec- 
tive contribution is of the order of (Av)-1. 

It is readily shown (see appendix I) that 


(12) sin 7 Av(t —r/Av) sinm Av(t — s/Av) 
_,  mAv(t—r/Av) x x Av(t —s/Av) 


if rss, and sono cross-terms will be found belonging to 
different sinc functions in the evaluation of the mean 
square electric field. But the rate at which energy is 
detected by a receiver is proportional to just this 
quantity. It therefore follows that the contribution 
due to any one sinc function term (or temporal mode, 
as we shall call it) is independent of the contribution 
due to any other mode. 

The integrated contribution to the mean square 
field due to one mode is 


o * gin? w Av(t —r/Ay) 
eer f [= Av(t —r/Av) 


A small receiver will draw from this temporal mode 
an amount of energy proportional to €,€ = But €, 
is just the electric field present at time LTA: 
and the statistical properties of the squared modulus 
of this quantity are described by equation (7). The 
amounts of energy Y absorbed by the receiver from 
the different temporal modes therefore obey a rela- 
tion analogous to (7), namely 


_ (13) (W2)=2(F), 
when classical physics applies. . 

The formula (13) has to be modified when quan- 
‘am physics applies. It is shown in Appendix II that 


dt =0 


aye e,e* Ay. 
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relation (13) is then replaced by ¥ =n hv, where n 
takes integral values, and obeys the relation 

(14) (n®)=(nj)+2(n}?. 

This formula can easily be modified to apply to the 
experimental apparatus of Hanspury Brown and 
Twiss. Consider two receivers S, and S, which are 
so close together as to be exposed to coherent radia- 
tion fields (in other words they jointly form a small 
receiver). Relation (14) for the joint receiver then 
states that 


( (7% + M)?) = (ny + ng) + 2(n, + ns )* 
or 
(ni) + (mg) +2 (mma) = (rm) + 2(m)? + 
+ (Mg) + 2 (me)? + 4 (m1) (Me), 
and in virtue of the formula (nj) = (n,) + 2(n,)? 
it then follows that 

(15) 

This is the basic result in the theory of the HANBURY 
Brown interferometer, and it is seen that it follows 
equally well from the classical formula and from the 
quantum formula. 

Only one generalisation remains to be made. If the 
electric fields E, and FE, at the receivers S, and S, 
are not perfectly correlated, so that 

FE, =TE,+ (1—|T |} Ey 
where F£, and £; are uncorrelated, then 
E,E* =" |£,|? + (1—|T|*)* EE 


(My Me) = Any) (n*) . 


, Say, 


and 

EES EES = 11" |£,|¢ + (4 —|T*) |£, |? | 2, |? 
16) 2 S(T ee ree 
The last term on the right hand side of (16) has the 


expectation value zero, since there is no phase rela- 
tion between EF, and £3. 


We write E,E* = U,, etc., and then find that 
(OU Ves (ICO) a OP Pray nC eae 

( U, )=( U?) =( U, )it follows that 
(U,U0,)=(1+ 19 12) ( U1) (Us }. 


Since 


|I'\? has been interpreted by Zernike [15] as the 
visibility 9 of fringes produced by the interference of 
light from S, and S,. In terms of n, and n, the for- 
mula reads 


(nyns)= (1 IT (7) (mm) | 


(17) = (1 + ¢) (m4) ( Me). 


* We have already seen that it makes no difference whe- 
ther classical or quantum theory is used to derive formula (15) 
for the coincidence rate in the arrival of photons at two 
small receivers exposed to coherent radiation. By analogy 
the classical result is here extended to include the case where 
the coherence is merely partial. 
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This completes the preliminary statistics. We 
finally note that 7 Av wavetrains pass in a time 
interval 7, and { n) T Av quanta are expected to be 
registered in this time. On the other hand, if there is a 
radiation flux of J photons/cm?/(c/s)/s, the expecta- 
tion value of the number of quanta detected in a 
time interval T by a receiver of equivalent area S and 
quantum efficiency xis x SJ T Av, where Av is the 
frequency bandwidth of the incident radiation. Com- 
parison shows that (n ) =x SJ. 


III. An Analysis of the Measurements. — The des- 
cription in this section relies heavily on the experi- 
mental results due to ForrESTER, GUDMUNDSEN and 
Jounson [5], who established the fact that quanta 
eject electrons from photo-cells with a time delay 
which is significantly less than 10-49 s. This means 
that the release may be considered to be practically 
instantaneous in the case of the present apparatus, 
whose amplifier has a bandwidth of the order of 
4108 els. 

We shall first consider a radiation field with a uni- 
form intensity J photons/cm?/(c/s)/s in the frequency 

Ay 
range (s mG, 
outside. The field is decomposed into temporal modes, 
a typical one having its maximum intensity at time 
r/Avy, and having appreciable intensity only in a 
short period around this instant. Suppose that the 
absorption of one photon causes the release of charge q. 
If the receiver at P absorbs n,(P) photons from the rth 
mode, the current released will be qn,(P) 8(t —r/Av). 

The current released during the interval from 
t = 0 until t = T is due to 7 Av temporal modes and 
equals 


(18) 


Ay) F : ' 
3 ) and with zero intensity 


5 Si 


in(t) = 9 n,(P) 8[¢ — (r/Av)] 


I IM 2 


where we have put N = 7 Ay. 

It is more convenient to express the current in 
terms of its frequency components, by means of the 
relation , : 

3) =_R +} eno een? | | | 
valid for the period (0, 7). The truth of the relation is 
easily shown on multiplying both sides by f(t + 7) 
and integrating over the interval (0, 7). This leads 
to an expression of f(t) in terms of its FouRIER 
series. 

In the complex notation 


MLS at (6) = ae » n,(P) i+ SS eae | 
i r=1 s=T 
See BOG ke Se 
= a > et “ oy ¢,(P) e—2mist/7 
n=] s=!1 
where 
(20) ¢(P) = > n.(P) e2tisr/N 
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The first term on the right hand side of (19) repre- 
sents the direct current component, which is filtered 
out. The alternating component is passed through a 
filter circuit with attenuation o, and phase shift ©, | 
at frequency s/Z. The current transmitted from P| 
becomes | 

Qt) i) =52 Do, ty 


$=1 \ 


j,(t) (t) —i(27st/T—e,) ; | 


with a similar expression for the current from Q. The 
integrated hee of the real currents is 


| 
y| 
Mes o 
he 
—= 
af 
_~ 
2) 
“~~ 
wx 
ox 
© 
ae 
m9 
a % 
Be) 
— 
wX 
= 
© 
— 
| 
eo ty 


s=1 

One expects to measure this signal. It will be confu- 
sed by random noise (shot noise) from the two recei- 
vers, whose effect we shall estimate shortly. We 
shall first find the value of &. 

Suppose that the cells have quantum efficiencies x 
and equivalent areas S, and S,. According to the — 
results found in § 2, 


Here e is the correlation between intensity fluc- 
tuations at P and Q. Wave-trains with different. 
suffixes are independent, so that 


(n,(P) n,(Q)) = (7,(P)) (7,1(Q) ) 
(25) = ny Ni 
when r=<r’. It follows that 
CPOs >| 


(7 CP ) een Tes cas eae de | 

(24) [MM(Q)) = ng = SQ J ; 
and (n,(P) n,(Q)) = (1 + 9) my ng | 

= (1+)? S,S, J? 


| 
> n,(P) n(Q) Fae | 
fae 


I 
" 
a | 


1 r=1 =| 


o - » ae ete | 
(26) | 


= en,n, I Av 


(1) Terms like j, (0) Jg () go out on integration from Oto T, 
as is readily seen by substitution from formula (21) for ie 
and from the analogous formula for To: 
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when s=<0, for in this case the second term in the 
square brackets vanishes, because 


2 2TTis/N ser 
(27) e2Tisr/N = eo | (1 __ e277 is) 
a 4 —e2tis/n 


It follows from (23) and (26) that 


Gy = 2¢q? NN, = | o Ay 
s=1 


(28) Sees G* S08 fot ops Ay: 
(Here fp = equivalent spectral width of the filter 
ercuit = F yc.) 
s=1 


The same result is found by Twiss and Hanpury 
Brown [14] who make their analysis in terms of 
infinite wavetrains of different frequencies. 

The formula for § can be generalised to include the 
case where J is not uniform over the bandwidth of 
the visual frequencies. Here the variation in intensity 
of the radiation field may be represented by a step 
function, with steps of typical width 8v. There is no 
correlation between the different frequency compo- 
nents, so that the effects of the different steps can be 
considered independently. We then find 
(29) CEG SS auio! pend? oY. 

Ay 

The effective duration of a temporal mode due to one 
step of width 8v is of order 1/8v. Our procedure 
seems valid if the apparatus cannot resolve so small 
an interval. This implies that Sv >> fp, or that the radia- 
tion field shall have no appreciable variations of 
intensity within a band of width fp, so that the step 
function shall be a good approximation. 

The sum over $v in the expression (29) for § can be 


replaced by an integral 


Ay 


J? dy, and the formula 


re-written with the aid of some auxiliary terms. 
Let Jia, be the maximum photon intensity within 
the band of light admitted. We define the equivalent 
bandwidth by 

‘| oe 
J dy 


(30) ghee 


max,/ Ay 


and the spectral concentration of the light by 


Acq i: J? dy | J? dy 
weiN oy 
(31) n, == > 4 2 =v 


|| ff ay] ae Avo 
Ay 


Phen, finally, 
(32) g = ey, 2 q? Doe fe NP ees Aveo: 
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It is to be expected that, in an experiment, the 
integrated product of the currents will differ from this 
value. The deviation is due to the fact that there is 
only a poor correlation between individual electrons 
released from the two cells owing to the low probability 
of absorption by either cell from any given temporal 
mode. The currents from the two cells therefore con- 
tain large uncorrelated components, and the integral 
for the expected product only begins to outweigh the 
random contributions after it has been building up 
for some time. We shall use the term « noise » for the 
product due to the uncorrelated currents, and shall 
find ‘Yb, the root mean square value of its integral. 
There are also terms due to the product of the random 
current from P, say, with the expected current from 
Q, but these are of smaller order and are ignored. 
We then have that 


t i he 4o ok ; : 
(33) w= (EP 0 Ro + ite iio at] 


e 


where it is now assumed that the fields producing 
the two currents are independent. 
Once again 


aq 


eo rake 
a4 fi, ih) + PO 10 4 
eo 0 
=F Daler) G@ + Fe) 4Q] - 


We are interested in the mean square of this expres- 
sion. In evaluating it we must remember that ave- 
rages can now be taken separately over quantities 
referring to the two receivers. The square of (34) 
contains terms like 


¢(P) t,(P) C*(Q) t%(Q) and %,(P) Ch(P) ¢,(Q) CQ) , 
whose expectation values are 
(PRC AP) PO). ON 
(ETP) al Cy Once Oe 

It will now be shown that 
(25) CP) CFD) a= | Gee Py aera yea 
and that 

(86) 202 (PYLE) Op whens Sea ce 


and 


* 


For WEP) se(b) 
N N 
aa N ai : (P) 5 (P) e2ti(rs—r's')/N 
pS r r 
r=1 7r'=1 
(37) 
Se ae ere et ar 
— ie 
Now (75.(P)ae(P)\ == mi(P)) (nih) = 254 
wasn, gGPeaae’ S Hence 
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hi ] eamir(s—s')/ po 
CAP) CEP) = [Pipe hy eee! ree mm Ss ae 
SAP NCAR) 2! (P)) — nz] (42) To =\/2 xg? (S, S,)* Jmax Mveq { Sof) 
\s= / 


N N 
i> a! t 

es N ~ ne e2tilrs— s \/N , 

a howd P 

r=a 74 

d ae 

The single sum has the value zero when ss’, 
equals 


N [( ni(P)) — nh 2 | =TAv|(ni(P))—n,], whens =s'. 


(38) (7) 


and 


The double sum is always zero, unless s = s’ = 0, as 
is easily verified by summing separately over r and 7’. 
In a similar way it can be shown that ( ¢,(P) ¢,/(P) ) = 9 
for all s and s‘, except s = s’ = 0. The proof is iden- 
tical, except that + s’now replaces —s‘ throughout. 
It follows that the expectation value is now zero 
unless s = —s’, or s = s’ =0 (since both s and s’ 
are non-negative). The components corresponding to 
s = 0 have, however, been eliminated from 7, and /, 
by the filter. The only surviving terms in the expres- 
sion for Tb” (equations (33) and (34) ) are 


Ye'= oe S of (¢,(P) C(P) ) ( &(Q) SQ) 


According to the results, given in § 2, 
(n2(P)) = n,(1 + 2 n,) 
for a small receiver ; we find that 


40) (n(P)) —n2 =n,(1+n,) =xS,J+x5,J), 


and 


TW = /2 xgX(S, S,)* J(L + xS, J)? (1 + x5, JP 


(41) 
- ‘ 1 5 
<= V2. xg2(S, S,)' E eg ped | 


(2 st) Jf UNS) 
s=1 

The extension used in the case of § can also be 
applied here to givea formula for Tb when the intensity 
is not uniform over the bandwidth of the radiation 
field. We simply add together the contributions from 
the different parts of the spectrum, finding, once 
again, that 


=, 
ps J 8y =| Jdv = Une Avo 
Ay a Av 


(7) Formula (38) can also be used to find the power spectrum 
of the output current j, (t). 


The factor in square brackets in (41) has a value 
near unity for low intensities, and has been omitted. 
The term © of 

sf 
spectral concentration of the filter circuits, defined by 


may be expressed in terms of the 


Y 4 ~ 4 
Dy a BS ace 
c= ==) 
(43) eee 
2 
o fou: 


od s 
s=1 


A comparison of (382) and (42) then shows that the 
ratio of the expected signal to the root mean square 
deviation in its measured value is 


rie fy 4 
Vi: 2249/2 (SS Te a) eee 
Vy 


to a sufficient degree of accuracy. The result, so far, 
holds for polarized light only. The extension to natu- 
ral, unpolarized light is readily made. Let the unpola- 
rized radiation, intensity J‘), have two equal, inde- 
pendent components with perpendicular directions 
of polarization, and intensities J! and J+. Then 

\ eae Le 

| and sl = yh 

in each frequency range. 

In this case the signal ¥ arises only from interac- 
tions between currents due to either the parallel or the 
perpendicular components. It follows that, in the 
present case, by formula (32), 

gg) = gl Ay gt 


SS ay) apt E (Jhex + Fx) 0 Meg 


(44) 9: 


(46) Sar = ug? S 54 fo r un aoe e Av, eq* 

As regards the noise, the number of photons pro- 
ducing it is simply doubled. Now we have neglected, 
in formula (42), the factor occuring in square brac- 
kets in (41). This is equivalent to neglecting the pos- 
sible multiple absorption of photons from any one 
mode and assuming, correctly, that the overwhel- 
ming majority of photo-electrons are produced singly. 
It is then easy toextend the formula for Yb to the case 
of unpolarized light, for there are now just twice as 
many quanta to eject the photo-electrons. Since Yb is 


linear in J,,,,, it follows that 
Vol — 9/2198 (S5Q)* (thax + Thax) A%eq (foT)* Vm 
‘47) et 

== 2'ngh (SS e)8 iE ay eG 
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Finally 


(48) yg) Ke 


Mss (SS) 


V/ max 
Ye 


The signal is largest when the incident radiation is 
most concentrated and y, = 1. On the other hand the 
noise is least important, for a given total bandwidth, 
when the electronic system has a small spectral con- 
centration. 

In a practical case one uses mirrors of, say, 2 square 
metres area and cells with 10 per cent quantum effi- 
ciency. The photon flux at 5 000 A from, say the star 
Sirius is J‘) = 2 x 10-* photons/cm2/(¢/s)/s. With 
an equivalent bandwidth f, = 108 c/s, spectral con- 
centrations 1/2 for both the light and the electronic 
system and one hour’s observation, we find 


(49) 


Giot )t Q: th. 


gt) sa) — 12 9:14. 


If there is perfect correlation between the fields 
at P and Q, so that p = 14, the signal : noise ratio is 
12:1. If the correlation factor | I'| is only 0.1, so that 
e = 0.01, the signal : noise ratio is 12 : 100. To obtain 
a ratio of 12: 1 the results need to be integrated for 
10 000 hours — far too long in practice. The appara- 
tus is therefore most useful for studying the variation 
of p for small distances of separation, where the corre- 
lation is high. 
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Appendix I. — Formula (12) is equivalent to the 
statement 


4 “sin x(@ —r) sin (2 —s) an 
ie esp GES) tees 
Now 
(50) sin x(x —r) =(—) *sinx(z — s) , 


if r and s are integers. 


Hence 
ff sin? x(z — 1) 
i Seay =| 
;? jo (tr) (%—S) 
is =e } sin? r(v—r) sin? x(x —7) ae 
i —S$) _« “£—r x—S 
a af \ sin? x(z—r) sin? x(%—S) Be 
Fo) Jie } “x—r L—s 
in virtue of (50), and so, finally, 
of ae \ \ sin? u sin? uw ee ie 
r—=-8. — 0 / U J 


inless r = s. 
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Appendix II. — We now establish the formula 
(707) = (7) pe Are ye 
for the general case of a source of natural radiation. 
Let a small receiver be exposed to linearly polarized 
radiation, having a uniform intensity in a small fre- 


Ay Ay : 

quency range (% ens a =) and coming from a 
source whichis stationary, in the sense of§ 2. The radia- 
tion field may be divided up into temporal modes. We 
consider the statistics of the number of quanta absorb- 
ed by the receiver in the various modes. Let a par- 
ticular mode be centred on the instant ty. The electric 
field at that instant may be expressed as a combina- 
tion of the fields due to different smaller frequency 
bands lying within the full range 


( Ay Ay 
Wi 9 » Yo t mi 


Let E, be a typical contribution, then 


if there are M smaller bands altogether. The intensity 
of the radiation at time fy is proportional to 
M 2 

a y 
oe 


s=1 


EE* = a 


which has the expectation value 
M 
(EE*) = ¥ (B, Ee), 
s=1 
since there are no phase relations between contribu- 
tions from different frequency ranges. The relation is 
equivalent to the statement 


M 


ety ( 1, ) 


about the number quanta in different modes. The 
square of the intensity of the radiation is 


Mi 2 M 2 
*)9 NS Wipe 
Noe \s=a 
M M M M 
= NI \) V * * 
— » yas >; had oe Ey, fie Ee, 
Sa aS) S21 ye 


which has the expectation value 
(BE)) =| 3 eet 
i 


—— 


| M M 
Bot Ne SS One epee ae 
+ 2 SOE Ee) (Bo) 
s=1 u=1 | 
su 


100 Pes 


an 
s=1 
| M \e M 
Ow NY Hees Si ( EE*) 
Sse 6 ie) Risto il ae s 
EE | s=1 


In terms of the number of quanta per mode absorb- 
ed by a small receiver, this equation becomes 


| M \ ; om ih M 
t Nae e yeh iN) SENSE) \2 
in ie he pares DIA 


\s=1 Sa 


Now n, must be a whole number. A typical small 


frequency band, say the s‘® one, makes only a small 


contribution to the radiation field, and so the proba- 
bility of n, exceeding unity is vanishingly small. 


But n= =n, when n, = 0 and when n, = 1 ; it follows 


that (n? ) =(n). On the other hand. 
(n ; Rive 
(n,) =O | MM | , and Say ie | =O | VI ve | 
7 
Mo (eet | 
gh hae 


and this term may be neglected in (51). We then find 
that. 


ax \ Lan \ 2 


2 Ny | N n \ + 2 SS h. 
Za tend 8 


=(n)p2(nj 


(52) 


It is well known that this relation also governs the 
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distribution of quanta in a radiation field in thermo- 
dynamic equilibrium. As has been shown here, the 
measurement of ( n? ) and of (2) by asmall receiver for 
light from an unresolved, stationary source will always 
lead to values satisfying relation (52) whether the source 
is in equilibrium or not. A measurement of this kind 
gives no additional information about the nature of 
the source. This result is in contradiction with a sug- 
gestion due to Fetiaerr [3, 4], but agrees with the 
experiments of Hanpury Brown and Twiss [ 9 ]. 
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Der Versuch, die Beugungserscheinung wellen- 
theoretisch zu erklaren, fihrten Youne@ zu der An- 
nahme, dass von jedem Punkt des beugenden Randes 
eine neue Kugelwelle ausgeht, die mit der urspriing- 
lichen Welle interferiert. Diese Hypothese der 
Beugungswelle erweist sich nun als ungemein frucht- 
bar, bei der Behandlung einer grossen Zahl von Beu- 
gungserscheinungen, wenn auch bei der exakten 
Losung der Beugungsprobleme die Aufspaltung in 
eine primiére und einer Beugungswelle sich nicht 
immer durchfihren lasst. Diese Aufspaltung gelingt 
jedoch in allen Fallen, wenn die vorliegende Aufgabe 


der Kircunorr’ schen Losung zuginglich ist. Denn 
das Kincnourr’sche Integral setzt sich seinerseits 
exakt aus einer einfallenden und einer Beugungs- 
welle zusammen. Aber auch in anderen Fallen zeigt 
sich die Annahme einer Beugungswelle vorteilhaft, 
weil auf diese Weise die Rechnung gegeniiber der 
exakten Losung wesentlich vereinfacht wird. Aus die- 
sem Grunde stellt der Verfasser die Beugungswelle in 
den Mittelpunkt seiner Betrachtungen und vermit- 
telt durch eine eingehende Diskussion der Grundlage 
der Beugungstheorie ein tieferes Verstindnis fiir die 
Zusammenhinge zwischen den verschiedenen Auf- 
fassungen. Auf diese Weise ist das in vielen Zeit- 
schriften und Biichern zerstreute Material hier gesam- 
melt und unter einem einheitlichen Gesichtspunkt 
dargestellt. Man muss den Verfasser fiir diese miihe- 
volle Arbeit besonderen Dank aussprechen. 
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Investigation into the theory of prediction of the appearance of colours 
and its bearing on the theory of colour vision 


E. G. T. Wasser, 
University College for Girls, Ein-Shams University, Cairo, Egypt. 


SumMARyY. — The evisting theory of prediction of the appearance of colours under different viewing conditions are reviewed : It appears 


that no unique set of three fundamental stimuli proved suitable jor deriving prediction formulae of general applicability on the 
basis of the trichromatic theory and the coefficient law. Nor was it possible to use these principles to derive a unique set of three 
fundamental stimuli by analysing the matrices of transformation which were obtained by applying the method of least squares to 
colour matches. The first observation led Mac ApAm fo doubt the adequacy of the trichromatic theory and to put forward a multi- 
receptor postulate, with the reservation that however numerous the receptors may be they should ultimately merge into three chan- 
nels to the brain. The implications of the second finding were not fully worked oul, bul some investigators took it, too, as evidence 
of the inadequacy of the trichromatic theory. 

In this paper an attempt is made to derive the form of the relationship between the tristimulus coefficients of corresponding 
colours from considerations of the Colour Equation, thereby evading the intricacies of the theory of colour vision. Supporting evi- 
dence of the linearity of the relationship is obtained from an analysis of adaptation data procured by binocular matching. The 
geometry of the transformations is studied, and some interpretations bearing on the concepts of vision are discussed, 


SOMMAIRE. — On passe en revue les théories existantes prévoyant l’aspect des couleurs sous différentes conditions d’observation. Il 


semble bien qu’un ensemble unique de trois excitations fondamentales ne convienne pas pour déduire une formule de prévision 
d’ application générale basée sur la théorie trichromatique et la loi des coefficients. Il n’est pas non plus possible de se servir de ces 
principes pour déduire un ensemble unique de trois excitations fondamentales en analysant les matrices de transformation obte- 
nues par l’application de la méthode des moindres carrés aux mélanges de couleurs. Le premier point fit douter Mac ApAm de 
la suffisance de la théorie trichromatique et Vamena a introduire un postulat basé sur UV existence de récepteurs multiples, avec la 
restriction que quel que soit le nombre de récepteurs ils devraient finalement se fondre en trois & leur entrée au cerveau. On n'a pas 
complétement étudié les conséquences qu’impliquent le second point, cependant quelques chercheurs Uadopterent comme preuve 
montrant UVinsuffisance de la théorie trichromatique. 

Dans cet article, on essaie de déterminer — en tenant compte de lV’ équation de couleur — le lien existant entre les coefficients 
trichromatiques des « couleurs correspondantes ». On évite ainsi les complexilés de la théorie de la vision colorée. On obtient la 
preuve du lien linéaire en analysant les données d’adapiation fournies par un systéme de mesure a vision binoculaire. On étudie 
les transformations géométriquement, et on discute de quelques-unes des interprétations portant sur les concepts de la vision. 


ZUSAMMENFASSUNG. — Ein Studium der bisherigen Theorie der Vorhersage des Farbenaussehens unter wechselneden Sehbedingungen 


zeigt, daB es kein einheitliches Tripel von Grundvalenzen zu geben scheint, das sich zur Ableitung einer allgemein anwendbaren 
Vorhersageformel aufgrund der trichromatischen Theorie und des Koeffizientensalzes eignet. Auch war die Anwendung dieser 
Prinzipien zur Ableitung eines einheitlichen Tripels von Grundvalenzen durch Auflosen der Transformationsmatrizen, die 
man mittels der Methode der kleinsten Quadrate bei Farbvergleichen erhdlt, unmdglich. Die erslere Beobachtung fithrle Mac 
ApvAm zu Zweifeln an der Brauchbarkeit der trichromatischen Theorie und zur Aufstellung einer Forderung von mehr Rezeptoren 
mit der Einschrdnkung, daB, so zahlreich die Rezeptoren auch sein mégen, sie schlieBlich uber drei Kandle im Gehirn einmiinden. 
Die Folgerungen der zweiten Entdeckung wurden nicht véllig ausgewertel, aber einige Forscher betrachteten sie ebenfalls als 
einen Hinweis auf das Versagen der trichromatischen Theorie. 

In der vorliegenden Arbeit wird versucht, die Form der Beziehung zwischen den trichromatischen Koeffizienten korrespondieren- 
der Farben aus Betrachtungen der Farbgleichungen abzuleiten, unter Vermeidung der Schwierigkeiten der Theorie des Farben- 
sehens. Eine Analyse von Adaptionsdaten, durch binokulare Farbenmischung gewonnen, unterstiitzte dies durch augenschein- 
liche Linearitét der Beziehung. Die Geometrie der Transformation wird untersucht, und einige Deutungen im Hinblick auf die 
Vorstellungen des Sehens werden erdrtert. 


I. Introduction. — To develop formulae for the 
prediction of the appearance of a colour it has been 
common practice to set out with some hypotheses con- 
cerning the mechanism of colour vision and the way 
the light receptors are affected by the adapting con- 
ditions. Most prediction formulae are based on the 
trichromatic theory of vision and the von Kries Coef- 
ficient Law; but prediction on this basis was not 
entirely successful. Some investigators attempted to 
deduce from the discrepancies evidence of the exis- 
tence of a multiplicity of receptors which obey the 
von Kries law. However, the fact that only three 
physical colours are necessary to match all colour 
sensations forces itself on the promoters of the multi- 
zeceptor theory : they suggest that the receptors merge 
into three channels to the brain. 

In this paper the current theories of prediction are 
sritically reviewed ; and an attempt is made to deve- 
iop an approach that is based on the less cuntrover- 
ial representation of a match by the colour equation 
ind the geometry of the transformation in place of 


hypotheses concerning the mechanism of colour vision. 
The implications on the theory of colour vision are 
then briefly discussed. 


Il. Methods of prediction based on the trichromatic 
theory and a discussion of the multi-receptor hypothe- 
sis. — A brief review of recent prediction formulae 
will elucidate the current trends : 

In 1952, Hetson, Jupp & Warren [1] derived a 
set of equations giving the changes in the appearance 
of colours when passing from daylight to incandescent 
filament illumination. The derivation of the formulae 
was based on the assumption that three independent 
photo-sensitive elements (or receptors) exist in the 
eye and respond to the red, the green and the violet 
regions of the spectrum respectively (The Younc- 
Hetmuoitz hypothesis). It was further postulated 
that the changes in the sensitivities of the receptors 
due to changing the state of adaptation in the eye can 
be taken into account by correcting the responses of 
the receptors by multiplying them by the factors that 
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would make the white colour always correspond to the 
excitation of the receptors in equal proportions. Let r,, 
£.) % be the tristimulus values of a colour sample S, 
illuminated by the illuminant C and measured in 
terms of the red, the green and the violet fundamental 
responses in the eye, and let r,, g, and ¢, be the corres- 
ponding coefficients of the colour when illuminated by 
the illuminant A. To make 7,, g,, %, (which refer to 
the colour of the sample to an eye adapted to illu- 
minant C not to illuminant A) correspond to the colour 
as it appears to an eye adapted to the latter illumi- 
nant, they should be multiplied by the factors k,, 


kg, k 

k.. a rclv a ’ ke == &cléa ’ 
where ra, Za, %4 and Tc, gc, % are respectively the 
tristimulus coefficients of the iluminants C and A mea- 
sured in terms of the fundamental response system 
of the eye. 

The theory was tested by a number of investigators 
using data obtained by different matching techniques, 
and employing various sets of values of the funda- 
mental responses. The main tests were : [41] HeELson, 
Jupp & Warren (loc. cit.) used the three fun- 
damental responses derived from the co-punctal 
points of colour defectives. The predicted values 
of colours under illuminant A were compared with 
memory matches. The agreement was poor in 
respect of the saturation of the colour samples and 
was not entirely satisfactory in respect of the hue. 
The authors suggested that the erroneous prediction 
of the saturation was probably due to the fundamen- 
tal system used, but they questioned whether it was 
worthwhile to make computations based on other sys- 
tems ‘in view of the approximate nature of the theory”. 

(2) Adopting the Wricutr-Tnhomson fundamental 
respose system [2], SanpERS [3] applied the theory 
to some data which the present writer [4] had obtai- 
ned in 1952 by means of a binocular matching method. 
The discrepancies between the predicted and the 
observed colours were ‘“‘ much larger than could be 
accounted for by the errors of observations ”’. 

(3) The third attempt to test the theory was made 
by Mac Apam [5] using data obtained by means of 
a local adaptation technique. Various combinations 
of the systems of fundamental responses were tried 
out, but ” large systematic discrepancies between the 
predicted and the observed colours showed up ”’. This 
led MacApam to inquire into the adequacy of the 
trichromatic theory and the von Krirs Law. Mac- 
ApAM examined the distribution of what he called 
‘the chromaticities of the discrepancies”, 24, Yq, 
defined as follows : 


wy = (XX G(X’ — XG) + (¥"— YG) =e(Z Ze) 
ya (YY + (yy te 
where X', Y', Z' are the tristimulus values of the 


observed matches and X¢, Y¢, Zé are the tristimulus 
values of the predicted matches. MacApam argued 


e where 


k, = TORAN ? 
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that if the discrepancies were due only to using the 
wrong values of the fundamental responses, the points 
which represent the chromaticities of the discrepan- 
cies on the chart should be the projective transforma- 
tions of the chromaticities of the original colours. As 
this was not the case he drew the conclusion that the 
visual process could not be explained on the basis of 
only three receptors in the eye. Some alternatives 
were examined. For instance, if there were more than 
three receptors in the eye and the prediction equations 
did not include corresponding terms, the chromaticities 
of the discrepancies would be randomly distributed 
on the CIE chart if the wrong values were taken for 
the first three fundamental resposes but certain trends 
should show up if the right values were used. In this 
latter case the chromaticities of the discrepancies 
would cluster round a point if there were four recep- 
tors, but would line up if the number of missing 
receptors was two. Three missing receptors would give 
a random distribution. Both the random distribution 
and the linear arrays appeared in the cases examined 
by MacApam. This led him to postulate that five or 
six receptors might be operating, although they 
might ultimately merge into three channels to the 
brain. 

It should however be noted that the interpreta- 
tions given of the distribution of the chromaticities 
of the discrepancies on the chart are not exhaustive. 
It is for instance evident that the inevitable errors 
of observations and rounding off would scatter the 
chromaticities of the discrepancies randomly if the 
trichromatic theory were correct and the right values 
of the fundamental respoises were used ; but in this 
case the discrepancies themselves would have been 
small and there would be no point in caleulating their 
chromaticities. On the other hand the clustering of the 
chromaticities of the discrepancies about a point, 
which MacApam interprets as indicative of a fourth 
receptor, may also be obtained if the prediction refe- 
rence system happens to be related to the correct 
system by the equations 


Ye=rY’, : 
Lie ie 


For in this case the chromaticities of the discrepan- 


cies cluster about the point 


Bay CL py) Ue ae) 
Also, if the relationship of the systems is of the form 
1 , 1— p)(1— 
ogame obec ue ae ety 
Yo == 1 ee 
Lue es 


the chromaticities of the discrepancies will array 
themselves close to the line given by the equation 


Xqy — Yq = constant. 
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The observed trends are not therefore conclusive 
evidence of the inadequacy of the trichromatic theory- 
A possible cause of the trends may be the use of a 
unique set of the fundamental responses for different 
pairs of states of adaptation (See § VI). 


Ill. The form of the relationship between the tris- 
timulus coefficients of corresponding colours derived 
from considerations of the colour equation. — The 
form of the relationship between the tristimulus coeffi- 
cients of corresponding colours may be seen from con- 
siderations of the colour equation without involving 
the number or the mode of operation of the colour 
receptors. 

Let us consider a colour C seen by a dark adap- 
ted eye and matched monocularly by the quantities 
r, g, 6 of the three instrumental primary colours R, 
G and B which are the red, the green and the blue 
colours respectively. The colour equation of C is 


(1) (C)y = r(R)p + g(G)p + 4(B)p. 


The valucs r, g, b are considered to be the coordinates 
of the colour C in the coordinate system R, G, B. If 
the state of adaptation of the eye is changed to another 
state L we shall still need the same quantities r, g, b to 
match the colour C, since a colour match remains a 
match when the adaptation of both the viewing and 
the measuring eyes are equally changed (cases of 
high-intensity adaptations are excepted). The colour 
equation representing the match in the new state of 
adaptation will therefore be 


(2) (G)ras (Py Pa(G);, +’ O(B)y: 


But if we view the colour C with an eye under the 
state of adaptation L and the matching mixture R, G, 
B with the other eye dark adapted, the match will be 
disturbed and we shall need different amounts 7’, g’, b’ 
of R, G, B respectively to re-establish the match. The 
colour equation representing the match will therefore 


be 


(3)  (C)p = 7'(R)p + 8'(G)p + 4'(B)p. 
Equation (2) represents the colours timulus which 
evokes the same sensation as G when the measuring 
primaries are viewed by the light-adapted eye, whereas 
equation (3) represents the colour stimulus that evokes 
the same sensation when the measuring primaries 
are viewed by the dark-adapted eye. The coordinates 
r, g, b andr’, g’, 6’ in the R, G, B system represent 
the “ corresponding colours ”,i. e., colours that evoke 
similar colour sensations if seen under two different 
states of adaptation. Thus, corresponding colours may 
be represented by the equation 


(4) r(R), + g(G), + (B)y =r'(R)p + g'(G)p + 5'(B)p- 


The red primary seen by the light-adapted eye and 
lenoted by (R);, does not of course match the same red 
orimary seen by the dark-adapted eye. The colour 
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(R')p which it matches can be represented in terms 
of the system R, G, B in the dark-adapted eye by the 
equation 


(5) (R’)p = o,(R)p + Bi(G)p + ¥1(B)p. 


Similarly, the colours G’ and B’ which match the 
green and the blue primaries seen in the light-adapted 
eye and denoted by (G);, and (B);, may respectively 
be represented by the equations 


(Cs Re a(R)p5 + Be(G)p + ¥2(B)p 
(B')p = ot3(R) 5 ai 23(G)p of y3(B) p. 


Equation (4) may thus be written 


(6) rR’ + eG’ + bB’=r'R + g'G + OB , 


where we have assumed that a colour in the light- 
adapted eye may be replaced by its match in the state 
of dark adaptation. Substituting for (R‘)p/, (G’)p and 
(B’')p we obtain 


rR+gG+0B=  r(aR + 6,6 + vB) 
+ g(%R + B.G + 2B) 
+ b(a;R + BsG + y3B) . 
Hence, 
r' = or + ag + agd 
g’ = Bur + Bog + B3d 
bf yal Yak =e Y3? 


which may be written 


ia Oy Ole Oe iy 
«| -|see |e]. 
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The coefficients of the colour C when seen by a light 
adapted eye are therefore linear functions of the coeffi- 
cients of the corresponding colour under the other 
state of adaptation. Since the CIE tristimulus values 
X, Y, Z of the colour can be obtained from the R, G, 
B coefficients using the linear transformation equations 
between the two systems, it follows that 


Xx! Xx 
| | =| 7] 
Li Z 


where Misa3 X 3 transformation matrix. 

The reader will at once notice that these equations 
are identical with formulae (7b) in BREWER’s paper [ 6 ] 
entitled ‘‘ Fundamental response functions and bino- 
cular colour matching ”. BrEewer’s theoretical treat- 
ment, however, proceeds from a different set of axioms 
involving the presence of three independent receptors 
in the eye which are affected by adaptation in accor- 
dance with the extended form of the coefficient law. 
Brewer’s approach will be taken up in a later sec- 
tion (§ V) ; we shall first seek empirical evidence of the 
linearity of the transformation. 


IV. Empirical evidence of the linearity of the trans- 
formation equations. — In recent years there has 
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accumulated ample evidence of the linearity of the 
transformation under varied experimental conditions. 


(1) The linearity of the transformation was tested 
using the binocular matching data which the author 
had obtained in 1952 [4]. The adapting illuminants 
were Standard Illuminants A and C. Illuminant A 
was obtained by tungsten lamps matching Standard 
IJluminant A as standardized at N. P. L., London; and 
the illuminant C was obtained by Cuance O. B.9 glass 
filters with illuminant A. MuNsELL samples were 
used as reference and matching samples. Interpola- 
tion and (to a limited extent) extrapolation were allow- 
ed between the discontinuous steps of hue, value 
and chroma of the samples of the MunsELt Book of 
Colour. Every match was the mean of four observa- 
tions, two by each eye. The reference samples, viewed 
with an eye adapted to illuminant A, and their matches 
as viewed with an eye adapted to the illuminant C are 
given in Table 1. The CIE tristimulus values of the 
matching samples were obtained from the O.S. A. 
Tables of the MunseL_i Samples. For details of the 
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apparatus and the observational procedure reference 
may be made to Vol. 2, N°3, of this Journal [ 8 }. 

The elements of the matrix of transformation were 
calculated by the method of least squares to make 
the sum of the squares of the residuals a minimum. 
The available thirty matches gave 


BOD —.148 19s 
M = | —.131 1 507 .068 |- 
—.197 .596 oe! 


The observed and computed matches are shown in 
figure 1. In general, the discrepancies are small and 
can be attributed to the following sources, in addi- 
tion to the errors of observations : 


(a) The CIE tristimulus values for the observed 
matches were calculated by interpolation from the 
tables. 

(b) The standard illuminant C which was used in 
the experiments was a combination of glass filters 
and illuminant A. 


TABLE 1 


Observed & predicted matches under illuminant C 
(White background) 
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Fig. 1. — Observed +- and predicted o matches (author’s data). 


(c) It is also possible that the colour perception of 
one eye was slighly affected by the state of adapta- 
tion of the other [ 4]. 

(2) The work of Burnuam, Evans and Newuatt|[ 9 ] 
on the ‘ prediction of colour appearances with diffe- 
rent adaptation illuminants”’ furnishes some _ evi- 
dence on linearity. In this experiment corresponding 
colours were determined by means of a binocular 
matching technique using three adapting illuminants 
A, G and G. The transformation matrices M,, M,,..., 
M, were calculated for the six combinations of the 
three illuminants. Each matrix was derived from 
twelve equations by the method of least squares. 
Acceptable predictions were reported in the cases 
(A, C), (C, A) and (G, C), but the criteria of acceptance 
are not given in the paper. 

(3) MacApam [5] mentioned that the evaluation 
of the transformation matrix from the data them- 
selves by means of the method of least squares gave 
better prediction than the use of Jupp’s theory with 
pre-determined values of the receptors, but the errors 
remained large and systematic. A possible source of 
discrepancy may however be sought in the local- 
adaptation technique itself, for the adaptation of one 
part of the retina to an illuminant may influence the 
solour perception of the part of the retina which is 
‘dapted to the other illuminant — especially since 
‘MiacApam used some coloured adaptations. Further- 


more, in MacApam’s technique the test and matching 
colours and the adapting fields were viewed alterna- 
tely by the eye. This is not quite the natural way of 
seeing colours, as colours are usually seen against 
their adapting background. 

To sum up, the analysis of the available binocular 
data seems to indicate that the tristimulus values of 
corresponding colours are related to each other by 
linear transformations. 


V. The Search for the Three Receptors. — Since 
the Younc-HEeLMHOLTz theory was proposed there 
has been an incessant endeavour to discover the three 
fundamental receptors by means of which the various 
phenomena of vision can be explained. An intensive 
search, using adaptation data, was initiated by 
Wricur [10] some years ago. Wright had hoped that 
“by measuring the colour changes induced on a 
given test stimulus by adapting the retina to various 
colours the results could be so analysed as to lead to 
the derivation of equations (for the three fundamer- 
tal resposes) in terms of the colorimeter stimuli R, G, 
B and, thence, ... the physiological response curves ”’. 
When this analysis failed, an alternative approach 
was attempted. In this the problem was inverted 
from a study of the effect of varying the adaptation 
while the test stimulus was unchanged, to an investi- 
gation of the effect of varying the test colour using a 
constant adaptation. The aim was to find three colour 
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stimuli in the red, green and blue regions that would 
not change in chromaticity when altering the state 
of adaptation of the eye, i, e., colours of invariant 
quality. Such colours should correspond to the sought 
receptors. 

By examining different monochromatic test colours 
with a chromatically adapted eye, Wricut discovered 
two regions in the spectrum where the test colours 
underwent minimum colour changes ; one region was 
near the red .70 w spectrum colour, the other was 
near the blue .46 vu. Wright accepted the first as the 
red fundamental and maintained that the blue fun- 
damental was not very far removed from the other. 
No monochromatic green colour was however found 
to remain invariant. Further assumptions were there- 
fore introduced to calculate the position of the green 
receptor ; but on the whole its determination was not 
satisfactory. 

A more generalized approach was adopted by 
Watters [11] who attempted to calculate the posi- 
tions of both the red and the green receptors from an 
analysis of the observations using an extended form 
of the coefficient law. The red stimulus was located 
close. to Wricut’s value, but the range of values 
obtained for the green stimulus did not contain the 
mean of Wricut’s determinations. 

The general theory of this approach was recently 
formulated by BREWER [6] whose method of analysis 
should have permitted the simultaneous evaluation 
on the three fundamentals. Starting from the asump- 
tion that the spectral response functions 7, g, 0 under 
on state of adaptation are related to the spectral res- 
ponse functions 7’, g’, b' under the other state of adap- 
tation by the equations 


g’ =k — ky 

where the coefficients k,, k,, k,, are specific constants 
of the particular pair of adaptation conditions, the 
conclusion was reached that the tristimulus coeffi- 
cients X, Y, Z and X’ Y’ Z’ of the corresponding 
colours are related by the equations 


XM axa per o¢ 
Wy, ff Tis VA 


where M-! = A-! K A, and 


Oo 
g° 


== 
r=ky, ? 


k, Ay, Ayg Ayg 
Ki ke and, A = | a1 Ag. a3 
hy, Q31 Age, Ags 


The matrix A is the transformation matrix that gives 
the response functions of the fundamentals in terms 
of the spectral reponse functions of the CIE Standard 
Observer. If the hypothesis of three independent 
receptors were true then the matrix A should be unique, 
and it can be determined from an analysis of adapta- 
tion data obtained under any pair of adaptation con- 
ditions. 

Analytical treatment of some matrices M obtained 
by Brewer [6] and Burnuam, -Evans_ and 
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Newua.t [9], as well as from MacApam’s data, gave 
no unique set of values of the fundamental responses. 
It is clearly futile to attribute the discordances to 
errors of observations, since in many cases the linear 
transformation gave acceptable predictions. A search 
must therefore be made for some loophole in the 
basic assumptions. 


VI. Geometrical Interpretations of the Linearity 
of the Transformation and its Implications on the 
Theory of Colour Vision. — The facts that emerge 
from the previous discussion of the theory and the 
practice of the prediction of the appearance of colours 
in a light adapted eye are. 


(1) That the tristimulus coefficients of correspon- 
ding colours are often related by homogeneous linear 
equations. 

(2) That the analysis of the adaptation data in 
accordance with the hypothesis of three independent 
receptors obeying the coefficient law did not lead to 
the determination of a unique set of fundamental 
stimuli, in spite of the fact that the linearity of trans- 
formation follows from such an assumption. 


Let us therefore consider the implications of the 
first finding. Geometrically, any linear transformation 
of coordinates leaves either one or three directions 
unaltered. The effect of the transformation on points 
lying on these directions is to displace the points along 
the lines radiating from the origin to their initial posi- 
tions. In terms of colour, these points represent colour 
stimuli that do not change their chromaticities in the 
new state of adaptation. To find the unaltered direc- 
tions of a transformation M let X, Y, Z be a point on 
one of these directions and let us denote the new posi- 
tion of the point by »X, AY, AZ, where 


AX = m4, = ipee & a M432 
AY = myX + MY + MZ 
WA = Tips -- Taek + M332 


Hence, 
100 Pes My2 My3 
Mei Mg: — A M3 = 0 
M34 M32 M33 — A 5 


which is a cubic in 4 and has either one or three real 
roots. Substituting each value of A into the transfor- 
mation equations and solving for X, Y, Z we obtain a 
set of values of the ratios X: Y:Z. Each set definies a 
direction that remains unaltered under the transfor- 
mation. 

The algebraic treatment of the matrix of transfor- 
mationis similar to BREWeER’s, but the basic reasoning, 
and consequently the interpretations of the results 
are essentially different. Brewer starts with the 
assumption of three independent receptors obeying 
the coefficient law in its extended form. The linearity 
of the transformations follow from this assumption. 
The values of A are the coefficients k,, k, and ky, and 
the sets of ratios X : Y: Z give the positions of the 
fundamental stimuli. If on the other hand the linearity 
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of the transformation is accepted as an empirical fact‘ 
or on the strength of the reasoning given in CSvItl ys 
the solution merely gives the colours which possess 
invariant quality when changing the state of adap- 
tation. 

The invariant colours corresponding to some pairs 
of adaptation conditions were computed from the 
transformation matrices M derived from the author’s 
own data [4] by the method of least squares. The 
matrices given by Burnnam, NewHatt and Evans te} 
were similarly analysed, with the exception of those 
matrices which did not give satisfactory predictions 
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of the colour appearances since the linearity of the 
transformation is a pre-requisite of the validity of the 
analysis. Table 2 gives the values of the unaltered 
chromaticities. The estimates of the fundamental res- 
ponses obtained by Wricut [10], Warters [11], 
Jupp [12] and Brewer [6] are included in the Table 
for comparison. The results are also plotted in figure 2. 
One of the invariant colours seems to fall in the red 
region and the other falls in the blue region. The third 
invariant colour, however, favours no particular 
region ; it appears everywhere in the chart, depending 
on the conditions of observation. 


CO-EFFICIENT OF Y 


CO-EFFICIENT OF X 


Fic. 2 a. — Invariant colours in the red and blue regions. 
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TABLE 2 


Invariant colours & fundamental stimuli 


Invariant colours : 
(1) (2) (3) 
2 y z x y z x y 
BREWER 2) Cla tAC rion erate cherie 741 248 —.025 . 754 .026 . 220 223 —.057 
BURNHAM, EvANS & NEWHALL : 
OS Cer Rn ee Comey IN OIRD 733 316 —.049 .958 — _ .053 095 .195 ss) 
GSAS Meiers ha ei nattertcene ote 732 316 —.047 1.018" —=7 6 098 ~19D .115 
A WES ope oes 8 Gin cern eoiente 912 —.069 2157 |e 3.64 — 1.33 214 -123 
WASSEF : 
1 ELANCE hie eee am Oto 700 .3/5 —.074 .924 — .664 . 740 149 041 
tee S i antanvan Aare HOD Ue 702 345 —.047 complex roots 
(Black backgd) 
Fundamental stimuli : 
IW RGEU Te staan cetay i atone Nh 4 Gr sleye Od .264 —.001 2.913 — 1.77 —1.14 .134 022 
(chromatic adapt.) 
WAMEIRS sya amet octets ete 749 .253 —.002 14.39 —12.47 —1.9 
(chormatic adapt.) 
SUID IIE yay apencen sy aitey nas Cis ae ke oe yce slant .747 .203 000 1.00 000 1.00 .16 000 
(dichromatic copunctal to to to 
points) 1.20 —.20 .18 
10 fall in the green region. This suggests that the third 
receptor in the eye is of anessentially different nature, 
shifting its maximum sensitivity with the state of 
adaptation. This would at once explain the failure o 
Wricur’s first approach to evaluate the fundamental 
05 stimuli from binocular data gathered under different 
states of adaptation. 
The present work, however, can not reveal the 
nature of the third receptor. Whether it is dependent 
9 on the responses of the other two receptors, or is of 
complex independent nature can only be decided by 
more experimental work appropriately designed to 
correlate the changes of the position of the third inva- 
~0-5 riant colour with the state of adaptation of the eye. 


Fia. 2 b. — Some locations of the third 
invariant colours. 


VII. Conclusions and Future Work. — The clus- 
tering of the invariant colours obtained under diffe- 
rent sets of adapting conditions in two regions on the 
chart, the findings of Wricur of the least colour 
changes in the same two regions, and the agreement 
of the results with Jupp’s location of the fundamen- 
tal from colour defective data lead to the conclusion 
that these regions may be associated with two receptor 
systems which are affected by adaptation according to 
the coefficient law. As regards the third invariant 
colour, it is noticeable that its place on the colour chart 
changes with the state of adaptation of the eye. It will 
also be remembered that there has always been mark- 
ed discordance between the various determinations 
of the third receptor process that was throught to 


The occasional presence of complex roots in the 
cases when linear transformations fit the data well 
calls for explanation. Individual differences between 
observers should also be introduced as factors in fur- 
ther experimental work to find out their bearing on 
the determination of the invariant colours. 

It is also important to examine the cases where the 
linear transformation does not adequately represent 
the changes of the appearance of colours. 
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Den wissenschaftlichen Beitrigen voran steht 
ein Aufsatz tiber Leben und Wirken von Orro EppEn- 
STEIN anlisslich der 80. Wiederkehr seines Geburt- 
stages. (Von H. Borcrnotp u. E. WaNDERSLEB). 

An rein optischen Beitrigen sind zu nennen : 


A. SONNEFELD, « Kritische Betrachtungen zu den 
Bedingungen iiber die Verzeichnungsfreiheit der 
Fernrohre » (S, 26-29), die sich mit dem umstrittenen 
Problem der Airy’ schen Tangentenbedingung ausei- 
nandersetzen. Der Aufsatz von H. RresenperG « Das 
Spiegelmikroskop und seine Anwendungen » (S. 30- 
78) beschreibt die praktisch wichtigsten Typen von 
Spiegelmikroskopen, insbesondere Abwandlungen des 
SCHWARZSCHILD’schen Zweispiegelmikroskops, und 
gibt an Hand einer Analyse ihrer spezifischen Bildfehler 
ilinweise fiir ihre sinnvolle Verwendung. Einer wellen- 
optischen Betrachtung tiber den Einfluss der Zen- 
tralabschattung auf die Intensitétsverteilung im 
Beugungsscheibchen und zweckmassige Beleuchtung 
folgen praktische Anwendungen im UV und IR, fir 
die ausschliesslich oder doch vorteilhaft Spiegelmi- 
kroskope in Frage kommen. Umfangreiches Litera- 
turverzeichnis (89 Zitate). 

Von H. Pontack sind zwei Beitrage enthalten. 

1) « Zur Theorie der absorptionsfreien achroma- 
tischen Lichtteilungsspiegel » (S. 79-86). Bericht tiber 
ein schon friiher von ibm verdffentlichtes Verfahren 
zur theoretischen Berechnung dielektrischer Inter- 
ferenzschichtsysteme fiir verlustfreie, farbsticharme 
Strahlenteiler, sowie Vergleichung mit  praktisch 
erzielten Ergebnissen. 


2) « Uber die reflexionsmindernde Wirkung diinner 


Metallschichten auf Glas » (S. 87-93). Uberblick tiber 
die Eigenschaften diinner metallischer Schichten auf 
Glas, sowie eigene Untersuchungen mit dem Ziel, 
Metalle ausfindig zu mache, die ein gutes Reflexions- 
vermégen besitzen, das ausserdem von der Wellen- 
lange moglichst unabhingig ist. M. Rremann « Uber 
die Verinderungen optischer Systeme durch die 
Verkittung von Linsen » (S. 216-228). 

Schon lange ist bekannt, dass Kittflachen die Qua- 
litt optischer Systeme oft stark beeintrachtigen. Die 
interferometrisch ausgemessenen Deformationen der 
Wellenflichen (zahlreiche Abbildungen im Text) wer- 
den erklart durch Verspannungen der Linsen unter 
Jem Druck der Fassung auf die Linsenrander, wobei 
“er Kitt zur Mitte der Kittflache ausweicht. 


E. Jxknicu, « Zur Theorie der Keildistanzmesser 
eater besonderer Beriicksichtigung eines neuen Distanz- 
vesserkeiles der VEB Jena » (S 140-215). 


Ernst Lav und WoxtrcanG Kru. Die Aquidensito- 
metrie, Grundlagen, Verfahren und Anwendungs- 
gebvete. 99 S. mit 106 Bild. Akademieverlag Berlin 
1956, Preis geb. DM. 19.50. 


Das Aquidensitenverfahren ist zunichst fiir die 
Auswertung von Interferenzstreifen entwickelt wor- 
den, hat sich dann aber auch bei der Beurteilung 
von Photogrammen mit kontinuierlich verlaufendem 
Intensitatsabfali bewihrt. Es stellt eine zweidimen- 
sionale Photoimetrie dar, wodurch erst die Photo- 
vratuioe ibren gesamten Informationsgehalt herge- 
ben. Fur die Gewinnung der Aquidensiten sind ver- 
schiedene Methoden entwickelt worden. Bei der ein- 
fachsten davon, der photographischen, wird das 
Negativ mit einem danach hergestellten Diapositiv 
zur Deckung gebracht. Dann ist als Folge des nicht 
linearen Verlaufes der Schwirzungskurve oder eines 
zwischen Negativ und Positiv verschiedenen Gam- 
mawertes die Schwarzung nicht tberall konstant. 
Eine photographische Kopie dieser Kombination 
zeigt die Kurven gleicher verbleibender Durchlassig- 
kert mit tiberraschender Schirfe. Gentigt sie fiir die 
Auswertung noch nicht, so kann man das Verfahren 
wiederholen und Aquidensiten zweiter und dritter 
Ordnung herstellen. Die Schwierigkeit, die darin 
liegt, Negativ und Positiv sehr genau zur Deckung zu 
bringen, lasst sich durch Ausnutzen des SABATIER- 
Effektes wahrend der Entwicklung des Negativs 
selbst umgehen. Auch die Schwiarzungsplastik lasst 
sich fiir die Herstellung von Aquidensiten heranzie- 
hen, wahrend das elektronen-optische Verfahren 
den Vorteil einer echten zweidimensionalen Photo- 
metrie aufgibt. Bei der Bedeutung, welche die Aquiden- 
sitometrie zu gewinnen verspricht, ist ihre Darstellung 
von der Hand der Begriinder des neuen Verfahrens 
sehr zu begriissen. 


G. Toratpo pi Francia. La Diffrazione della Luce. 
Edizioni Scientifiche Einaudi, 1958, 530 pages. 


L’auteur fait une mise au point des résultats récem- 
ment acquis dans le domaine de la diffraction. I 
donne tout d’abord l’outillage mathématique indis- 
pensable (en particulier la transformation de Fou- 
RIER) et rappelle les lois classiques de la propagation 
des ondes. 

Dans un large tour d’horizon sur l’ensemble des 
phénoménes de diffraction que l’on rencontre en 
optique, l’auteur reprend l’étude des problémes anciens 
en les présentant sous une forme moderne et insiste 
sur les résultats obtenus dans ces dernieres décennies. 

Titres des chapitres : 

1te partie : Introduction mathématique. Les ondes 
lumineuses et leur propagation. 

2° partie: Les réseaux infinis a traits rectilignes. 
Les réseaux réels et leurs applications. Les dia- 
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phragmes non périodiques : a) unidimensionnels, 
b) bidimensionnels. 

3° partie : Diffraction et formation des images. Les 
images des objets lumineux. Images dobjets éclairés 
en lumiére cohérente. 

4° partie: Réseaux généralisés. Critéres de Vinté- 
gration approximative ; principe d’ HuycEns-FRESNEL. 
Application du principe d’Huycens-FRESNEL. Quel- 
ques paradoxes optiques. 

5° partie: Géométrie des ondes aberrantes. Ondes 


BIBLIOGRAPHIE 


[ Opr. Acra 


aberrantes et diffraction. L’interférométre a réseaux. 


Nous regrettons pour notre part qu’une place 
plus large n’ait pas été faite a la théorie de la forma- 


tion des images optiques, théorie qui touche peut-étre 


plus lingénieur spécialisé que le physicien. 


Cet ouvrage est incontestablement un document 


de grande valeur ; il nous semble devoir trouver 
une large audience et faire partie de la bibliothéque 
de fond du physicien. 

A.M. 
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